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Tre] 


Com um número de 1993, 

três números referentes a 1994, 

a Técnica conseguiu 

atingir 

aquilo a que se propôs: 

quatro números por ano. 

Razoavelmente espaçados, 

originais científicos de interesse, 

um noticiário diverso e rico sobre a Escola, 

uma publicidade atraente, fonte ainda de informação. 

Tanto bastou para que novo original 

afluisse, 

os convites fossem aceites, os artigos convidados aparecessem, 

que as comunicações científicas se tenham afirmado 

e diversificado. 

Com estes factos, 

motivos de esperança, 

se inicia 1995. 

E inicia-se já 

com quatro notas breves, 

tais como as pedimos relativamente 

à Ciência, Tecnologia e Desenvolvimento. 

Sentimos e sabemos interesse 

nesta modalidade, nova, de comparticipação: 

Notas sintéticas e rigorosas, 

permitindo aferir o que é novo, 

comparar, medir avanços. 

Serve a Técnica a Escola e a Comunidade Científica, 

generosamente ajudada pelo Conselho Directivo do I.8.T. 

que deste modo também ajuda a AEIST 

na realização 

de uma obra que perdurará: 

A Técnica — Revista de Engenharia, 

integralmente devotada 

à Investigação e Desenvolvimento a 
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Prestar serviços associados à garantia de qualidade e segurança das 
obras, para além do ensaio de novas tecnologias e sua aplicação em 
engenharia civil, são algumas das actividades desenvolvidas pelo 
LNEC — Laboratório Nacional de Engenharia Civil. Às novas tecno- 
logias são transferidas para os utilizadores através de acções 
de formação e divulgação, bem como de normalização e 
homologação. 


OQ patrimônio natural, o património construído, produtividade e acti- 
vidades de apoio à indústria da construção, são algumas das linhas 
que estruturam o actual plano de investigação do LNEC. As actividades 
de I&D prendem-se com o desenvolvimento e difusão de conhe- 
cimentos no espectro alargado da engenharia civil, em particular nos 
domínios das obras públicas, habitação, urbanização, recursos hidricos 
e ambiente, para além do aperfeiçoamento de produtos e processos 
para a indústria da construção 


Apoiar a elaboração de normas e regulamentos portugueses na area 
da engenharia civil, é outra das atribuições do LNEC, destacando-se a 
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contribuição que tem sido dada pelo Laboratório para o 
aperfeiçoamento da normalização no contexto comunitário e 
internacional, 


Uma parte significativa das actividades é desenvolvida por meio de 
cooperação com outras entidades nacionais ou estrangeiras, atraves 
de conventos ou protocolos. 


[em assumido um particular relevo a cooperação com a União 
Europeia e os serviços prestados ao abrigo de convénios celebrados 
com instituições brasileiras e dos Países Africanos de Lingua Oficial 
Portuguesa, 


A aposta na cooperação é, de resto, uma forma de o LNEC evitar uma 
expansão excessiva dos meios humanos ao seu serviço, sem 
deixar de investir cada vez mais na respectiva formação. 


Sob a responsabiliddae do LNEC, encontra-se o Laboratório de 
Engenharia Civil de Macau, que constitui um instrumento 
priviligiado para a cooperação com este território e com a China. 


The progress in physics cer- 
tainly will depende to a large 
extent on the progress of nonli- 
near mathematics, of methods 
of solving nonlinear equations. 
W. Heisenberg” 


RESUMO 


Faz-se uma introdução à ciên- 
cia dos sistema complexos 
através de exemplos simples, 
escolhidos de modo a realçar 
alguns aspectos das relações 
entre o não-linear, o colectivo, 
o caos, a previsibilidade e a 
calculabilidade. 


A CIÊNCIA DOS 
SISTEMAS COMPLEXOS 


Rui Dilão 
Departamento de Fisica do IST 


1. INTRODUÇÃO 


Quando se junta um grande número de sistemas, as propriedades macros- 
cópicas ou colectivas do sistema composto não estão, em geral, relaciona- 
das com as propriedades dos seus constituintes individuais. Neste caso, o 
sistema composto é um sistema complexo. 

As ciências como a física, a biologia, a química, a economia, a história e a 
medicina estudam sistemas complexos: têm que lidar com as catástrofes, 
com a turbulência, com as doenças, com as revoluções, com a evolução na- 
tural, com a extinção dos dinossaurios, com a evolução da bolsa e do uni- 
verso. Com o desenvolvimento de máquinas de cálculo poderosas que che- 
gam a realizar 10'' operações de vírgula flutuante por segundo, é possível 
cálcular estes sistemas. 


A grande dificuldade no estudo das propriedades dos sistemas complexos é 
de que, em geral, os modelos matemáticos associados conduzem à determi- 
nação de soluções de equações não-lineares, sendo difícil a aferição de so- 
luções numéricas com os resultados experimentais e com as soluções anali- 
ticas (quando existem). Como a experiência tem mostrado, aparecem 
dificuldades relativas à fraca previsibilidade que muitos algoritmos fornecem, 
tendo-se criado a necessidade de desenvolver técnicas específicas para a 
análise de sistemas não-lineares. A teoria dos sistemas dinâmicos tenta 
cumprir este programa. 

O objectivo da teoria dos sistemas dinâmicos é estudar as propriedades 
qualitativas das soluções de equações, estabelecendo novos métodos de 
analise, observação e descrição dos sistemas reais. A teoria qualitativa dos 
sistemas dinâmicos ocupa-se das propriedades topológicas das soluções de 
equações. A teoria métrica dos sistema sdinâmicos ou teoria ergódica es- 
tuda o problema da possibilidade de determinar valores médios de grande- 
zas associadas à descrição dos sistemas. 


Enquanto os métodos quantitativos estão por excelência associados a pro- 
blemas lineares, a teoria dos sistemas dinâmicos desenvolve técnicas quali- 
tativas de analise global de soluções de equações com o objectivo de deter- 
minar e classificar os seus tipos genéricos. Esta abordagem permite prever 
propriedades dinâmicas, como sejam o aparecimento ou desaparecimento 
de singularidades, a emergência de propriedades estocásticas, etc.. Recen- 
temente, por análise directa de séries temporais obtidas experimentalmente, 
estas técnicas qualitativas têm vindo a ser utilizadas na previsão da evolu- 
ção de sistemas cujas equações dinâmicas não existem ou estão mal defini- 
das!. 

O mundo dos sistemas complexos é o mundo das transições de fase, das 
mudanças bruscas de comportamento ou das bifurcações, da emergência 
de ordem em processos caóticos, da ocorrência de processos violentos 
como a explosão de uma supernova ou a extinção de uma espécie. Como 


4 DINDAL 


1 


56 dNW | 


estes processos estão invariavelmente associados a modelos matemáticos 
não-lineares, esta exposição está organizada tendo como ponto de partida 
as relações entre o linear e o não-linear, tentando mostrar a diversidade de 
fenómenos que podem surgir nos sistemas complexos. 


Nos exemplos que se seguem, iremos ver como se relacionam alguns des- 
tes problemas, realçando as relações entre o não-linear, o colectivo, o caos, 
a previsibilidade e a calculabilidade. 


2. SISTEMAS LINEARES VERSUS SISTEMAS NÃO-LINEARES 


Nos sistemas de equações lineares, sejam elas equações às diferenças, 
equações diferenciais ordinárias ou equações às derivadas parciais, é sem- 
pre possível determinar soluções que obedecem a condições iniciais e con- 
dições fronteira. Mais ainda, através de um conjunto de soluções particula- 
res, é sempre possivel gerar novas soluções através de combinações 
lineares das soluções dadas. 

Esta ultima propriedade designa-se por princípio da sobreposição, técnica 
utilizada na construção de soluções gerais de equações, e tem aplicações 
importantes na mecânica, no electromagnetismo e na mecânica quântica. 
Por isso é muito frequente associar esta propriedade à própria estrutura e 
ao âmbito das teorias em questão. 


Dos vários domínios da física, o electromagnetismo e a mecânica quântica 
são particularmente sensíveis a esta discussão. Como é conhecido, as 
equações que determinam os campos e as funções de onda são lineares, e 
o principio da sobreposição aparece sempre como uma ferramenta para a 
descrição dos fenómenos que são tratados nestes domínios. 


Por exemplo, Werner Heisenberg, referia-se à mecânica quântica da se- 
guinte maneira”: 

Linearity in quantum theory has a very deep, almost philosophical reason 
and is not just connected with some approximation. In quantum theory we 
do not deal with facts but with possibilities: the square of the wave fune- 
tion describes the probability, and the superposition of the wave function 
the possibility of adding two solutions to get a new solution, is absolutely 
essential for the whole foundation of quantum theory. Therefore it defini- 
tely would be wrong to say that the linear character of quantum theory is 
approximate in the same sense as the linearity of Maxwell's equations is 
approximate. 


Por outro lado, em relação ao electromagnetismo e à electrodinâmica quân- 
tica a relação entre o linear e o não-linear é completamente diferente. Max 
Born e Leopold Infeld escreveram”: 


In all these cases there is sufficient evidence that the present theory (for- 
mulated by Dirac's wave equation) holds as long as the wave-lengths (of 
the Maxwell or the de Broglie waves) are long compared with the “radius 
of the electron” e /mc”, but breaks down for a field containing shorter 
waves. The non-appearence of Planck's constant in this expression for 
the radius indicates that in first place the electromagnetic laws are to be 
modified; the quantum laws may then be adapted to the new field equa- 
tons. 

e mais tarde, Barbashov e Clhernikov comentavam':: 


Since the appearence of the paper of Born and Infeld on a nonlinear elec- 
trodynamics of the free field it has become clear that the nonlinear theory le- 
ads to qualitatively new and much richer physical concepts than the linear 
theory. 


3. EXEMPLOS DE SISTEMAS NÃO-LINEARES COM COMPORTAMENTO 
COMPLEXO. EMERGÊNCIA DE PROPRIEDADES COLECTIVAS. 


3.1. A equação de Burgers: formação e propagação de singularidades 


A equação de Burgers descreve a evolução de um sistema de partículas 
sem interacções mútuas. É um o exemplo de como se podem formar estru- 
turas complexas (singularidades) quando se passa da descrição do movi- 
mento de uma partícula para a descrição (colectiva) do movimento de um 
conjunto de partículas. 

Como é bem conhecido, uma partícula isolada tem um movimento rectili- 
neo e uniforme, e a sua lei de movimento é x(t)=vyi+x, em que x, é 
a posição da partícula no instante t=-0 e "» é a sua velocidade. Vamos 
então considerar um conjunto infinito (contínuo) de partículas isoladas e 
sem interacções mútuas. No instante t=0, a distribuição espacial de velo- 
cidade e densidade é, respectivamente, volx)=v(2,t=0) e 
Polx)=p(x,t=0), figura 1. 
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Fig. 1: Distribuições iniciais de velocidade (vo (x) e densidade (Po(x)) de um conjunto 
de partículas. 


Como o movimento de cada partícula obedece à equação de Newton, x =0), 
a velocidade das partículas num instante arbitrário t é 


dv O dd dr A 


->— + — =— +yvr—=0 
de dt doa dt ú dx 


A equação de Burgers para o movimento colectivo das partículas é 


dt Ox (1) 


Como se mostra facilmente”'*!, a solução da equação de Burgers, obede- 
cendo às condições iniciais Volx) e Po(x) , é 


ví Xi F= Vo (x te tv(x,t)) 
PolXo) 


Plxt)= 
Ler (2) 


em que x, é a coordenada Lagrangeana, o = x(t=0), 

Na figura 2 estão representadas as funções velocidade e densidade para 
t=1,0 e t=1,5, calculadas? a partir de (2) e das distribuições iniciais re- 
presentadas figura 1. 
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Fig. 2: Evolução temporal de uma distribuição de partículas, ao fim dos tempos ! = |O e 
t=1.5 calculada a partir das distribuições de velocidade e densidade da figura 1. 


No instante += 1,0 eno ponto x=«, a derivada de v(x.!) em ordem a x é 
infinita e para +>1,0, existe uma região do espaço em que a velocidade 
das partículas pode tomar dois valores. Do mesmo modo, através do gráfico 
da densidade de partículas conclui-se que em 1=1,0 formou-se uma singu- 


laridade na função Plx.t): plx..t=1)=<>, Da análise exacta das soluções 
da equação de Burgers decorre que existem apenas soluções bem definidas 
para tempos 


Se em geral um sistema de equações às derivadas parciais não-lineares 
evolui para soluções que podem tomar mais do que um valor no espaço, 
surge o problema de saber como interpretar as soluções quando estas apa- 
recem associadas a sistemas reais. Existem duas propostas bem diferentes. 
Para alguns autores, baseados em resultados experimentais e na necessi- 
dade de previsibilidade de sistemas reais, as soluções singulares são apro- 
ximadas por soluções construídas através da regra da alavanca de Maxwell, 
figura 3. Neste caso, diz-se que a solução da equação desenvolve uma 


Fig. 3: Formação de uma frente de choque na equação de Burgers. A frente de onda de 
choque é determinada pela igualdade das áreas a) e b): Regra da alavanca Maxwell. 


onda de choque ou frente de onda de choque, a que está associada uma 
velocidade de propagação característica. A demonstração desta regra empi- 


rica deve-se a Hopf”. Contudo, do ponto de vista computacional e para sis- 
temas em que existe o controlo analítico do processo de formação deste 
tipo de singularidades, os resultados numéricos são diferentes dos resulta- 
dos teóricos. 


Outros autores preferem considerar que sistemas com um comportamento 
deste tipo, os vários valores admissíveis em cada ponto do espaço podem 
de facto ocorrer e, não existindo critério de escolha, os resultados experi- 
mentais são interpretados como sendo de natureza estocástica. 


De facto, ambas as interpretações são possíveis. Se se pretende seguir 
uma frente de onda de choque, a primeira interpretação é adequada. Con- 
tudo, quando se mede um campo ou uma velocidade no interior de uma sin- 
gularidade, os resultados experimentais vão a favor da segunda interpreta- 
ção. À primeira interpretação explica bem a formação de cogumelos e de 
jactos em explosões de estrelas”. A segunda interpretação dá conta, por 
exemplo, da observação de oscilações aleatórias do campo magnético da 
magneto-esfera, na região de choque do vento solar. 


Deste exemplo conclui-se que as soluções de equações às derivadas parci- 
ais não-lineares (tecnicamente equações quasi-lineares) podem gerar fenó- 
menos complexos como sejam a formação de singularidades e o apareci- 
mento de aleatoriedade durante a evolução temporal dos sistemas, 
tornando-os difíceis de prever. Quando se analisa o comportamento das so- 
luções de equações às derivadas parciais não-lineares e sobre as quais não 
existe controlo analítico, é difícil a calibração dos resultados numéricos 
quando comparados com observações reais”, 

Face a este panorama, é corrente afirmar que é possível existirem fenóme- 
nos naturais cuja possibilidade de previsão consiste na observação. 


3.2. À equação das ondas 


Uma pequena perturbação num meio continuo propaga-se de acordo com a 
equação 


o" 
ot” Ox” (3) 


== Co 


No entanto, esta equação é apenas uma aproximação à lei de propagação. 
De facto, quando se analisa o problema das vibrações de uma corda chega- 
se à equação não-linear 


To Mo 2 
Com a aproximação habitual das pequenas amplitudes, /( | +[58) ] =, 


obtem-se a equação das ondas (3), apenas válida para o estudo da propa- 
gação de pequenas perturbações. 

Na figura 4 está representada a evolução temporal de uma corda presa nas 
extremidades, calculada através da equação linear (3) e da equação não-li- 
near (4). Como se avalia facilmente por análise das figuras, o espectro da 
solução da equação linear mantem-se invariante ao longo do tempo, a me- 
nos de variações de amplitude, enquanto que para a solução da equação 
não-linear aparecem novas frequências. Levanta-se então o problema de 
saber como é que a energia se espalhou pelos restantes modos próprios. 
As simulações apresentadas na figura 4 não foram prolongadas no tempo, 
pois para +>(),155 a solução numérica da equação não-linear não con- 


verge, embora a técnica utilizada na obtenção das soluções apresentadas 
na figura tenha sido baseada na resolução formal da solução geral”. 


Equação linear Equação não-linear 


Fig. 4: Evolução de uma onda progressiva na equação de ondas linear (3) e não-linear (4). 


No caso de cordas relativistas (c = 1) a lei de propagação de uma onda 
obedece à equação de evolução” 


A aa apo a: 
dx) dr dt) dx (5) 


o” oO d dx diox 


Na figura 5 está representada a evolução temporal de duas perturbações na 
corda infinita mostrando a formação dinâmica de singularidades. Enquanto que 


Fig. 5: a) Formação de singularidades na interacção de duas perturbações numa corda in- 
finita relativista. b) Superficie integral correspondente à região do espaço em que se for- 
mam singularidades numa solução particular de equação (5). 


noutras equações não-lineares como a Korteweg-deVries a interacção de 
dois pacotes de onda do tipo solitão não gera nenhum tipo de estrutura na 
zona de interacção. No caso da equação (5), a solução qà(x.t) pode tomar 
vários valores nos mesmos pontos x e t, figura 5a). Isto é, se o observável 
d(x,t) pode tomar vários valores num ponto do espaço e não existe ne- 
nhum critério de escolha da boa solução, é de esperar que o resultado de 
uma observação seja uma sequência de valores aleatórios. 


As soluções representadas na figura 5a) foram obtidas analíticamente atra- 
vés da parametrização das soluções'” e não foi possível encontrar um mé- 
todo numérico que permitisse resolver as singularidades. O que é interes- 
sante reter neste exemplo é de que para pequenas amplitudes a interacção 
entre as duas frentes de onda não apresenta qualquer problema técnico, 
não se gerando singularidades, enquanto para grandes amplitudes apare- 
cem e desaparecem singularidades, todas elas parametrizada pelo tempo. 
Na figura 5b) estão representados dois cortes da superficie integral de uma 
solução da equação (5), mostrando o mecanismo de formação, evolução e 
desaparecimento das singularidades. 


3.3. Sistemas químicos. Formação de Padrões 
À reacção química autocatalítica 


k, avo ks 
B+X>' Ff 2X 4+Y =3º 3X (6) 


é descrita pelas equações cinéticas 


E rB(c-r)-H(C-Yf 
dt l 
Di) 

dt | 

X()=C-Y(1) 


(7) 


em que X, Ye B são as concentrações das espécies químicas e 4, e 
são as constantes das reacções. A constante C = X(r)+Y(t) é invariante ao 
longo do tempo, estando relacionada com a lei da conservação da massa. 
Com, k,=k,=1,0 e as concentrações iniciais X(0)=0,5 e Y(0)=0,0, é 
possível mostrar que as concentrações podem evoluir para dois estados de 
equilíbrio genéricos (X=C,Y=0,8=0) e (X=0,Y=C,8B2>0).0O primeiro 
estado é atingido se B<0,617.Se B>0,617, podemos ter uma infinidade 
de soluções de equilibrio. Na figura 6, representamos a evolução das con- 
centrações no espaço de fases (X.65B). 


Se esta reacção é efectuada num reactor agitado, garantindo uma homoge- 
nização dos vários componentes químicos, qualquer um dos estados de 
equilibrio do sistema de equações (7) é atingido nas condições indicadas. 
No entanto, se a reacção se dá num meio espacialmente extendido, como 
acontece por exemplo nos tecidos biológicos, as várias espécies químicas 
difundem-se, passando a existir competição entre reacção e difusão. Neste 
caso, as equações de evolução das espécies químicas são: 


OX ; Px OX 
— =-kBX+kRYX +D, : o 
dt dci | dE dy | 
oY x oy dy 
=kBX-kY X “+ D,| — 
do ú (E Í dy | 
dB “(B q9ºB 
% E | o ay | 


(8) 
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Fig. 6: Espaço de fases da equação cinêtica (7). 


em que D,, 2, e D, são os coeficientes de difusão no meio das três es- 
pécies químicas. Sistemas destes tipo, que têm soluções cinéticas elemen- 
tares, podem gerar fenómenos bastante complexos. Por exemplo, esco- 
lhendo D,=0,01, D, =0,0001 e D, = 0,00001, a concentração espacial da 
espécie química X evolui conforme se mostra na figura 7. 


t=10 t=50. 


Fig. 7: Evolução temporal do sistema de reacção-difusão (8), para condições iniciais alea- 
tórias e condições fronteira de fluxo nulo. Integração numérica com o método de Euler 
numa rede de 150x150. Ao fim do tempo 1=50 o sistema atingiu o estado de equilibrio. 


Assim, a reacção química (6) gerou estrutura espacial através de uma inte- 
racção difusiva, embora as leis cinéticas sejam elementares. Este tipo de 
acoplamento é um exemplo de como se podem formar padrões em sistemas 
bioquímicos, levando muitos autores a pensar que este tipo de fenómenos 
está na base da morfogénese"''?. Na figura 8 estão representadas duas dis- 


tribuições espaciais obtidas em reacções autocatalíticas '*. 


3.4. Caos determinista 


Quando se observa a evolução (temporal) de um sistema e se selecciona 
uma variável, por exemplo x, como representativa do seu estado, a hipótese 
de trabalho mais simples e a mais utilizada nas aplicações é de que a evolu- 
ção da média temporal de x(!) nos dá uma medida aproximada de como o 
estado do sistema está a evoluir. À teoria ergódica estabelece as condições 
de existência deste valor médio. Contudo, existem sistemas reais em que es- 
tes valores médios podem estar mal definidos, não assumindo nenhum valor 
particular. Neste caso estamos na presença de sistemas aleatórios, 


Fig. 8: Distribuição espacial de concentração obtidas em dois processos autocatalíticos di- 
ferentes, 


em geral, associados à formação de singularidades nas soluções das equa- 
ções dinâmicas, como é o caso dos exemplos 3.1 e 3.2. Nos sistemas com 
caos determinista, estes valores médios existem sempre mas a dinâmica do 
sistema é aleatória ou caótica. 


Para determinar se um sistema tem caos determinista existe uma metodolo- 
gia simples que consiste em avaliar como é que um sistema se comporta 
face a pequenas variações nas condições iniciais. Se, para pequenas varia- 
ções nas condições iniciais ao fim de um tempo t o estado do sistema 
afasta-se exponencialmente no tempo, na métrica do espaço de fases, en- 
tão o sistema tem sensibilidade às condições iniciais. Em geral, um sistema 
com sensibilidade as condições iniciais é caótico. 

O sistema caótico mais simples é o modelo conhecido por transformação de 
pasteleiro ou ferradura de Smale'*. Vejamos então como se constroi este 
modelo. 


Suponha-se um quadrado de lado 1 ao qual se aplicam as transformações 
(figura 9): 


1) Duplica-se a coordenada horizontal e divide-se por 2 a coordenada verti- 
cal. 


2) Sutura-se verticalmente em x=1 e empilham-se os dois rectângulos. 


Fig. 9: Protótipo de um sistema caótico: construção da transformação de Smale ou de pas- 
teleiro. 


Considerando que o quadrado inicial representa o conjunto de todos os va- 
lores possíveis das variáveis dinâmicas de um sistema no instante t=0, no 
instante t=1 o sistema evoluiu para um novo estado em que alguns dos pon- 
tos representativos do estado do sistema têm uma ordem relativa que foi al- 
terada de acordo com as transformações 1) e 2). 


Vejamos então o que acontece quando se pinta uma figura no quadrado. 
Como se vê na figura 10, os pontos coloridos vão mudando de posição ao 
longo do tempo e a figura coerente inicial desfez-se completamente. Muitos 
autores designam a transformação de Smale por transformação de paste- 
leiro (por razões obvias...). 
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Fig. 10: Mistura de condições iniciais e propriedades aleatórias da transformação de 
Smale. 


Ora, é possível mostrar que para esta transformação determinista, dois pon- 
tos muito próximos no quadrado inicial afastam-se exponencialmente com o 
tempo, e assim a transformação de Smale tem sensibilidade em relação as 
condições iniciais. Por outro lado, este sistema é equivalente ao processo 
estocástico do jogo da cara-ou-coroa, protótipo de sistemas com proprieda- 
des estocásticas. 


Todo o sistema dinâmico com esta propriedade é um sistema dinâmico caó- 
tico! 


Para encontrar sistemas que tenham estas propriedades basta então 
construir transformações com a capacidade de “dobrarem” ao longo do 
tempo o conjunto de estados iniciais. Ora, isto pode ser conseguido atra- 
vés de transformações de um intervalo, sobrejectivas e não invertíveis. 
Por exemplo, a transformação quadrática do intervalo [0,1], 
Co=dx(l-x |), figura 11a), tem esta propriedade. Na figura 11b) está 
representado a evolução temporal de um ponto genérico do intervalo 
[0,1], concluindo-se que a evolução temporal da transformação é aleato- 
ria, embora o sistema dinâmico seja determinista, isto é, um ponto do do- 
mínio transforma-se univocamente num ponto do contradomiínio. De facto, 
é possível mostrar que este sistema é caótico. 


a) b) 


X41 
ii 


X4 t 


Fig. 11: Gráfico da função quadrática do intervalo [0,1] e série temporal caótica. 


A transformação da figura 11 pode ser mergulhada na família de transforma- 
ções +,, =4ux,(I-x,), dependentes do parâmetro u, e May'?, num con- 
texto de dinâmica de populações, descobriu que ao variar o parâmetro 
u el0,1) podiam existir transições bruscas entre comportamentos assintóti- 
cos regulares periódicos & comportamentos caóticos (ou ergódicos). Na fi- 
gura 12 apresentamos o conjunto de estados assintóticos acessíveis da 


transformação x,, =4ux,(I-x,), em função do parâmetro 4. 


Fig. 12: Diagrama de bifurcações em função do parâmetro 4 para a função quadrática 
do intervalo. X,y = 44 x,(1—-x,) 

Até agora mostramos como surgem comportamentos caóticos em modelos 
matemáticos não-lineares. Seria interessante que eles aparecessem em sis- 
temas reais. De facto, num simples circuito eléctrico com uma junção p-n de 
capacidade variável foi possível observar o mesmo comportamento qualita- 
tivo que o encontrado na figura 12'%. Na figura 13 está representado o dia- 
grama de bifurcações a dinâmica caótica do sistema em função de uma ten- 
são exterior. As figuras 12a) e 12b) foram obtidas a duas temperaturas 
diferentes. 


Fig. 13: Diagrama de bifurcações obtido experimentalmente num circuito com uma junção 
p-n com comportamento não-linear. 


Outro exemplo de sistema caótico é o atractor estranho de Lorenz'”, associado 
à evolução temporal da temperatura no interior de um fluído, mantido entre 
duas placas paralelas a temperaturas diferentes. A equação de Lorenz é 


x=0(y—x) 
VY=Ix-y-xz 


2z=-bz+xy 
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em que x é proporcional à função de corrente e y e z são duas funções re- 
lacionadas com a temperatura. Os parâmetros re O são proporcionais, 
respectivamente, aos números de Reynolds e Pranditl e b é um factor de 
forma. Na figura 14 está representada a órbita no espaço de fases (x.).7) 
de um ponto inicial genérico. Neste exemplo, as órbitas de fase não conver- 
gem para nenhum estado estacionário, percorrendo erraticamente uma re- 
gião limitada do espaço de fases. 


Fig. 14: Atractor estranho de Lorenz no espaço de fases (x,Y.2), para »=28, 
o=I0eb=8/. 


Quando se varia o parâmetro r encontramos o mesmo tipo de estrutura que 
as encontradas nos mapas do intervalo. Na figura 15, representamos o dia- 
grama de bifurcações do atractor de Lorenz em função de r. Comparando 
este diagrama de bifurcações com o da figura 12, mostra-se a analogia entre 
o comportamento das soluções da equação de Lorenz com as soluções da 
equação discreta do intervalo, construída a partir da transformação de Smale. 
Conclui-se assim que o sistema de Lorenz tem comportamento caótico. 


Pensa-se que o comportamento encontrado no atractor Lorenz é comum a 
muitos sistemas não-lineares com pelo menos três variáveis de estado. 
A associação de processos caóticos a dinâmicas não-lineares deterministas 
abre um novo capítulo na ciência moderna, possibilitando a utilização de 
modelos matemáticos no controlo e previsão da evolução temporal destes 
sistemas. Em particular, é hoje um assunto de intensa investigação o de- 
senvolvimento de métodos de controlo não-linear para sistemas caóticos!”, 


4. CALCULABILIDADE E IMPREVISIBILIDADE 


Para terminar, vamos introduzir um sistema caótico simples mas cujas pro- 
priedades computacionais contradizem todos os resultados analíticos. 


seja a transformação do intervalo 


x =8(4,)=2h,|-| (9) 


. no 
20.0 SO ES r 100,0 


Fig. 15: Diagrama de bifurcações em função do parâmetro r para o atractor de Lorenz. 


com x, e[-1,1]. Ora, como o declive da função g(x) é +2,0s 2" pontos 
fixos das iteradas g"(x) são instáveis. Por outro lado, é possível mostrar 
que este sistema é caótico e tem sensibilidade às condições iniciais, o que 
implica que as iteradas de um ponto genérico do intervalo [-1,1] percorrem 
densamente esse intervalo. 


No entanto, quando se itera, por exemplo, o ponto *9 =0,97397, ao fim de 
53 iteradas a sua imagem cai no ponto fixo instável x=1, figura 16. Ora, 
isto contradiz todos os resultados analíticos e a própria afirmação de que o 
sistema é caótico. 


RES a 


Fig. 16: Orbita do ponto x, = 0,97397, por iteração da função do intervalo (9). 


Qual a razão desta discrepância entre resultados analíticos e observação 
numérica? A resposta é simples e está associada à representação digital 
dos números nos computadores. Na figura 17 está representado o gráfico 
da função continua (9), assim como o seu gráfico na representação interna 
de um computador a 4 bits. Devido à precisão finita dos computadores, pois 
estes trabalham apenas com números racionais, os pontos fixos instáveis da 
função (9) são transformados em pontos fixos estáveis e o sistema com- 
porta-se numericamente como um sistema regular, apenas com um longo 
transiente. 
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Fig. 17: Gráfico da função g(x) e da sua aproximação na representação interna de um 
computador. 


Terminamos com este exemplo, para chamar a atenção para o facto de que 
a previsão numérica de sistemas não-lineares pode fornecer soluções bas- 
tante diferentes das encontradas tanto nos modelos matemáticos como na 
realidade da experiência. Assim, os critérios de validação dos modelos nu- 
méricos tem de ser confrontados, por um lado, com os resultados das previ- 
sões teóricas (consistência), por outro, com os resultados experimentais (ca- 


libração). 
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HEINRICH HERTZ 


M. de Abreu Faro 


Original Submetido em: 4 de Janeiro de 1995 


Em 20 de Outubro de 1994 a Aca- 
demia das Ciências de Lisboa reali- 
zou uma sessão expressamente de- 
dicada à comemoração da vida e 
obra de Hertz, isto por ocasião do 
primeiro centenário da sua morte. 
Nessa cerimônia foram apresenta- 
das duas comunicações, uma pela 
Professora Lídia Salgueiro, outra 
pelo autor desta breve nota. Algu- 
mas passagens foram retiradas, na 
integra, da nossa comunicação na 
Academia que aliás é bastante mais 
extensa 


Heinrich Hertz nasceu em 22 de Fe- 
vereiro de 1857 e faleceu em 1 de 
Janeiro de 1894. 

À carreira científica de Hertz radica- 
se na Escola de Berlim, liderada por 
Helmholtz. Ai se dedicou na década 
de 70 especial atenção aos princi- 
pios fundamentais da Electrodiná- 
mica, ai se confrontavam as Teorias 
do Continente com a Teoria de 
Maxwell. 


Nesse espírito, Helmholtz propôs 
um prémio no âmbito da Academia 
Prussiana das Ciências. Que, em 
última analise, se investigasse ex- 
perimentalmente a existência de 
qualquer relação entre as forças 
(campos) e a polarização dieléctrica 
dos materiais isolantes. 

Para tanto, convidou, em 1879, 
Hertz que declinou o convite por du- 
pla razão: ser um trabalho difícil, 
vago, imprevisível que não se inte- 
grava num doutoramento que o 
queria mais simples e breve; não 
dispor de meios de medida que lhe 
permitissem evidenciar, com segu- 
rança, quanto se pretendia o que, 
de facto, tomava a questão ainda 
mais difícil. 

Declinou o convite mas a ideia 
ficou. 


Em 1886 obteve uma cadeira na 
Technische Hochschule de Karls- 
ruhe. 

Uma demonstração pedagógica 
com um jogo de bobinas planas, 
enroladas em dupla espiral — espli- 
rais de Riess ou Knochenhauer — 
revelou-lhe, indiciou-lhe, com sur- 
presa, possibilidades experimentais 
desde que recorresse, como recor- 
reu, ao fenómeno da ressonância. 
Isso permitiu-lhe dispor de uma 
sensibilidade insuspeitada 


Assim se iniciou uma série impor- 
tante de experiências que se 
iniciaram com o estabelecimento 
experimental de oscilações electro- 
magnéticas de frequências nunca 
antes e sistematicamente conse- 
guidas, da ordem de 10º Hz 
(100 MHz), o que foi publicado sob 
o título “Sobre as oscilações eléc- 
tricas muito rapidas”. 

No âmbito destas experiências des- 
cobriu experimentalmente o efeito 
fotoeléctrico, mais tarde explicado e 
teorizado por Einstein. 

O objectivo proposto por Helmholtz 
foi atingido, mas os acontecimentos 
precipitavam-se e para Hertz come- 
caram a configurar-se aspectos in- 
teressantes, genuinamente caracte- 
rísticos das teorias de Maxwell. 


Um deles era o facto do espaço li- 
vre de matéria se comportar como 
um material dieléctrico. 


A corrente de deslocamento no vá- 
cuo criaria um campo magnético. 
Era um facto que se tinha que acei- 
tar mas uma dificuldade subsistia, 


Hertz não conseguia demonstrar se- 
paradamente que, no ar, campos 
eléctricos variáveis criassem cam- 
pos magnéticos e que, também, 
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campos magnéticos variáveis crias- 
sem campos eléctricos. 

Mas .... e, para nós, é o pensamento 
mais fecundo de Hertz, pressentia 
que podia demonstrar simultane- 
mante os dois efeitos se conse- 
quisse provar a existência de ondas 
electromagnéticas propagando-se no 
ar com velocidade finita 

Assim se encontrou com o que teo- 
ricamente Maxwell previa em 1864 
quando, pela primeira vez, anunciou 
as equações gerais do campo elec- 
tromagnético e, simultaneamente, 
propôs uma teoria electromagnética 
da luz. 


Com esta motivação e após delica- 
das experiências atinge este último 
objectivo que excedia em alcance e 
profundidade o que em 1879 Helm- 
holtz tinha proposto. 


Surge assim em 1888 a sua comu- 
nicação fundamental “Sobre a velo- 
cidade finita de propagação de ac- 
ções electromagnéticas”. 

Ainda, em 1888, publica “Sobre as 
ondas electrodinâmicas no ar e sua 
reflexão”. 

Culmina os seus trabalhos com um 
trabalho teórico notável, diriamos 
que é o primeiro trabalho teórico 
que se fez sobre uma antena de 
emissão, foi publicado sob o titulo 
“As forças das oscilações eléctricas 
tratadas de acordo com a Teoria de 
Maxwell”, foi publicado em 1889, aí 
se exploram as potencialidades do 
“vector de Hertz”, por forças en- 
tenda-se campos. 


Em 1890 Hertz publicou dois traba- 
lhos teóricos notaveis: 


— Sobre as equações fundamen- 
tais da Electrodinâmica dos cor- 
pos em repouso. 

— Sobre as equações fundamen- 
tais da Electrodinâmica dos cor- 
pos em movimento. 


São trabalhos que, recorrendo à 
simplicidade do essencial, se de- 
senvolvem com uma lógica impecá- 


* Edmond Bauer - L'E'lectromagnetisme - Hier et Anjourd'Hui, Albin Michel, Paris, 
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vel, decorrente das hipóteses assu- 
midas. 


Para Hertz, a Teoria de Maxwell era 
o conjunto das Equações de Max- 
well. 


Essas equações traduziriam então 
leis experimentais que podiam ser 
testadas. 

As cargas eléctricas e as correntes 
eléctricas sao tomadas como sim- 
bolos matemáticos e representam, 
respectivamente, pontos limite de 
convergência ou divergência das li- 
nhas de força do campo eléctrico e 
eixos dos vórtices das linhas de 
força do campo magnético. 


Na exposição de Hertz não se re- 
corre a qualquer hipótese mecani- 
cista. 


Heriz entendia e salientava que es- 
peculações sobre a natureza da 
matéria impediam o progresso, 
eram questões ilegitimas. 

Coerentemente, Hertz, não discutiu 
a natureza da electricidade nem dis- 
cutia a essência do éter. 


Esta atitude havia de conduzir a di- 
ficuldades. 


Como sabiamente salienta Edmond 
Bauer”, a crise que se motivou com 
a Teoria de Maxwell, que foi desen- 
cadeada por Hertz e se continuou, 
noutra atitude, com Lorentz, revela 
que a Ciência, se quizer progredir, 
e do mundo real nos facultar con- 
cepções positivas e consistentes, 
terá, e quando isso for necessário e 
convier, passar do ponto de vista do 
fenomenológico a hipóteses de es- 
trutura e inversamente. 

Como nos diz M. Tonnelat**, a teo- 
ria electromagnética da luz, de Ma- 
well, nada nos diz nem garante so- 
bre o arrastamento parcial do éter 
da Óptica, teorizado por Fresnel. 
Como salienta Max Bom, Heriz foi o 
primeiro a utilizar a hipótese da 
completa convecção do éter pelos 


dieléctricos em movimento. 


“* M.A. Tonnelat - Histoire du Principe de Relativité, Flamarion, Paris, 1971. 


Hertz adoptou, assim, a hipótese 
formulada por Stokes para os fenó- 
menos da Óptica. 

Admitiu a completa convecção do 
éter pelos dieléctricos em movi- 
mento e isto no que respeita ao 
Electromagnetismo, 

Estava plenamente consciente da 
arbitrariedade que cometia. 


A teoria de Fresnel e a experiencia 
de Fizeau, e isso sabia-o Hertz, 
constituiam-se em indicação para 
que assim não fosse. 


A teoria de Hertz revelou fenomeno- 
logia nova, qualitativamente era tal- 
vez aceitável, quantitativamente 
não, como ficou demonstrado pelas 
experiências de Rôntgen (1885) e 
Eichenwald (1903), relativamente ao 
deslocamento lento de um isolante 
num campo eléctrico e, bem assim, 
pelas experiências de Wilson 
(1904), relativamente ao desloca- 
mento no seio de um campo mag- 
nético. 

A teoria de Hertz era simples e ex- 
plicava correctamente os fenóme- 
nos de indução em condutores em 
movimento. 


O seu último trabalho é de natureza 
teórica: “Os princípios da Mecância 
apresentados numa nova forma”. 

É uma obra póstuma, prefaciada 
por Helmholtz, que faleceu no 
mesmo ano de Hertz. 


Hertz não deu especial atenção a 
eventuais aplicações das ondas 
electromagnéticas. 

Hertz foi despertado por Helmholtz 
no ano da morte de Maxwell, 1879. 
Marconi, foi despertado por Righi no 
ano da morte de Hertz, 1894. 

Uma notícia necrológica extensa e 
descritiva da obra de Hertz, escrita 
por Righi, surpreendeu Marconi, de 
férias nos Alpes Italianos. Tinha 20 
anos. Nascera em 25 de Abril de 
1874. 


Volta a casa, à Villa Griffone, perto 
de Bolonha. 


Ai repete Hertz e vai além 
Parte para Inglaterra 


Em 2 de Junho de 1896 requer uma 
patente para a Telegrafia sem Fios, 


o que acontecia pela primeira vez 
no Mundo. 


Em Dezembro de 1901 consegue 
uma ligação radiotelegráfica da Eu- 
ropa à América. 

Em 1924, após intensas pesquisas 
experimentais a bordo do seu barco 
Elettra, telefonou, via rádio, de In- 
glaterra para a Austrália 

As ligações hertzianas directas, por 
difracção ou via ionosfera, eram 
uma realidade que marcou profun- 
damente o Século XX. 


Na realidade foi uma tríada notável: 


Maxwell - Hertz - Marconi 

Foi desta sequência que resultou a 
fama maior para Hertz 

Mas Hertz não foi só isso. 

Corajoso, não se submeteu ao pre- 


A 


conceito, a magestade do que ja 
era. 


luminou as suas experiências com 
o poder da sua observação e crítica 
científica 

Não experimentava ao acaso: pre- 
via, antevia, desistia, por vezes, pe- 
rante a indicação de sinais físicos 
tão fracos que não podia detectar. 
Assim e por isso foi declinado em 
1879 o convite de Helmholtz. 


Buscava e filiava-se nos princípios 
primeiros, ansiava por uma simplici- 
dade suficiente. 
Conseguiu experimentalmente o 
que Faraday já intuira e tentara. 
Disso sabia ele. 


Confirmou e contribuiu poderosa- 
mente para a difusão da Teoria de 
Maxwell. 


Consciente de que esta teoria não 
esclarecia nem tinha conseguido 
ainda aquilo que, de facto, a com- 
pletaria, preocupou-se com a Elec- 
trodinâmica dos Corpos em Movi- 
mento. 

Ai definiu uma linha de investigação 
onde avultou Lorentz e culminaria 
com Einstein: desinventor do éter e 
criador da Teoria da Relatividade. 

A nossa profunda admiração e res- 
peito para o físico ilustre, lógico, con- 
sistente e corajoso que foi e soube 
ser Heinrich Rudolf Heriz, erros que 
teve, foram, quem sabe, virtudes de 
inteligência, lógica e carácter. 
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PRÉMIO DA BOA | 
ESPERANÇA 1993 ATRIBUÍDO 
A PROFESSOR DO IST 


Dos dois contemplados pelo Prémio 
da Boa Esperança, um deles é o 
Professor Amaral Fortes do Insituto 
Superior Técnico. 

Este prêmio visa encorajar a activi- 
dade cientifica e tecnológica exer- 
cida em Portugal por cidadãos na- 
cionais ou estrangeiros e o 
aproveitamento dos seus resulta- 
dos para o bem-estar dos Portu- 
queses e de toda Humanidade. 


O trabalho de investigação e de- 
senvolvimento distinguido é um ar- 
tigo que foi publicado na revista 
Interface Science intitulado “Stabi- 
lity and Dissociation of Quadruple 
Junctions in Polycrystals”, 


Trata-se de uma análise termodi- 
nâmica sobre a evolução de jun- 
ções entre quatro interfaces que 
se intersectam numa linha. Esta 
configuração é geralmente instá- 
vel, dissociando-se em duas jun- 
ções de três interfaces (junções 
triplas), mas pode, em certos ca- 
sos, permanecer estável durante 
algum tempo até atingir uma contfi- 
guração em que a estabilidade não 
é mais possível. Estas junções for- 
mam-se transitoriamente na evolu- 
ção de estruturas celulares, por 
exemplo em policristais e em es- 
pumas. 


PLASMAS 
COLOIDAIS 


Prof. David P. Resendes 


Departamento de Fisica, .8,T. 
Coordenador Europeu da Rede “Plasmas Coloidais” 
Original Submetido em: 24 de Janeiro de 1995 


Plasmas coloidais constituem uma 
subarea emergente da física dos 
plasmas. O estudo destes sistemas, 
teve origem na fisica espacial com 
o estudo de vários aspectos da di- 
nâmica de grãos carregados numa 
tentativa de compreender alguma 
fenomenologia complexa ocorrendo 
em plasmas astrofísicos. Contudo, 
um grande numero de observações 
experimentais recentes revelaram a 
presença de grãos carregados em 
atmosferas terrestres e extra-terres- 
tres bem como nalguns equipamen- 
tos laboratoriais e industriais, Na 
astrofísica os grãos têm um papel 
importante na física, quimica e na 
dinâmica de nuvens interestrelares, 


no processo de formação de estre- 
las e planetas bem como na expan- 
são do invólucro de novas e super- 
novas. No sistema solar, plasmas 
coloidais ocorrem nos cometas e na 
magnetosfera de todos os planetas 
exteriores (Júpiter, Saturno, Urano e 
Neptuno). A interacção de anéis 
planetários com a magnetosfera 
(produção plásmica, absorção e 
transporte) bem como a ionosfera 
(injecção de massa e correntes) são 
fenômenos intensamente estuda- 
dos. A nível terrestre, tais grãos 
carregados são importantes na fí- 
sica e evolução da mesosfera, na 
interacção de aerossois na alta at- 
mosfera, na propagação de ondas 


electromagnéticas e na formação de 
relâmpagos. Investigação laborato- 
rial básica inclui efeitos observados 
em reactores de fusão, cristalização 
de suspensões coloidais, transições 
de fase e fenomenologia crítica. Em 
fins de Janeiro reuniram-se no Insti- 
tuto Superior Técnico em Lisboa os 
membros responsáveis da Rede de 
Investigação Europeia intitulada 
“Plasmas Coloidais” para definir as 
questões mais importantes e uma 
estratégia de colaboração concer- 
tada a nível europeu para os próxi- 
mos dois anos. À rede é constituída 
por seis laboratórios para além do 
IST, nomeadamente o Rutherford 
Appleton Laboratory (RAL), as Uni- 
versidades de Oxford, Napoli e 
Tromso, e Bochum, e o Max Planck 
Institute. Do IST, para além do au- 
tor esteve presente o Prof. Tito 
Mendonça do Dept. de Física. 
Neste artigo, serão muito sumária- 
mente abordados os avanços mais 
recentes e os problemas mais pre- 
mentes surgidos nesta reunião 
acerca desta nova área da física 
dos plasmas. 

Um plasma coloidal pode definir-se 
de um modo pouco preciso como 
um plasma com uma fase dispersa 
de objectos sólidos (grãos ou poei- 
ras). Tipicamente os grãos variam 
desde uma pequena fracção de mi- 
cron podendo ir até um centímetro 
de tamanho e com cargas eléctricas 
da ordem dos 10º a 10º cargas 
electrónicas. Plasmas coloidais dife- 
rem doutros plasmas em pelo me- 
nos três aspectos essenciais: 1) As 
cargas nos grãos não possuem um 
valor constante mas são determina- 
das pelas propriedades das próprias 
partículas e pelas propriedades do 
plasma envolvente. As cargas são 
portanto variáveis independentes do 
sistema. 2) Os grãos encontram-se 
envolvidos por partículas do plasma 
de modo que fortes efeitos não line- 
ares são importantes na vizinhança 
dum grão. Daí que grãos são sítios 
de fenómenos fortemente não linea- 
res e bastante mais estáveis que 
solitões, cavitões, etc. 3) Os efeitos 
colectivos produzidos pelos grãos 
nos plasmas possuem um tempo de 


vida longo o que significa que po- 
dem ser fontes eficientes de movi- 
mentos colectivos excitados em 
plasmas. 

Em astrofísica o efeito de grãos in- 
terestrelares na propagação de on- 
das de choque através de nuvens 
interestrelares continua a ser uma 
área activa de investigação. O 
amortecimento de flutuações de alta 
frequência, os efeitos na química in- 
terestrelar via interacções gás-grão 
e efeitos electromagnéticos media- 
dos pela inércia dos grãos através 
do campo magnético constituem al- 
guns dos fenómenos de interesse. 
Utiliza-se uma abordagem MHD 
multi-fluidos em que os grãos são 
tratados como diferentes fluidos in- 
teractivos (consoante a sua carga). 
Reacções químicas são incorpora- 
das no código. O efeito dos grãos é 
suavizar a estrutura da onda de 
choque originando temperaturas 
mais baixas. Um dos objectivos fu- 
turos é alargar o âmbito do trabalho 
para se tratar regiões gasosas de 
maior densidade para se poder des- 
crever choques de acreção durante 
as primeiras etapas de formação de 
estrelas e desse modo obter sinais 
identificativos destes eventos invisi- 
veis directamente por obscuração 
dos grãos interestrelares. 


Nas aplicações ao sistema solar 
destaca-se a interacção entre anéis 
planetários, em particular a compo- 
nente fina de grãos e o campo 
magnestosférico, plasma e ionos- 
fera. Neste contexto, há uma varie- 
dade rica de física a explorar: a fi- 
sica do processo de carga dos 
grãos, focagem electrostática, bem 
como processos plásmicos como a 
velocidade critica de ionização, 
“pick-up” de electrões e iões pelo 
campo magnético corrotativo, isotro- 
pização do plasma e termalização 
entre os vários componentes. Para 
além destes processos microscópi- 
cos, estão a ser abordados também 
processos colectivos como o arras- 
tamento do tubo de fluxo associado 
a correntes de grãos. A importância 
deste efeito nas magnetosferas dos 
planetas exteriores será abordada. 
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A nível terrestre, prossegue-se a 
uma análise de plasmas coloidais 
na presença duma componente 
neutra significativa. Processos fisi- 
cos como fricção, gravidade, coli- 
sões entre iões e componentes 
neutros, bem como recombinação a 
superfície de grãos e o efeito fotoe- 
léctrico estão a ser sistematica- 
mente estudados. De elevado inte- 
resse é o processo de carga de 
grãos, gotículas de água e flocos de 
neve - estudos aplicáveis às áreas 
de difusão de ondas em nuvens e o 
relampejar atmosférico. No contexto 
industrial (e laboratorial) geram-se 
grãos ou poeiras em descargas de 
rádio frequência utilizadas no pro- 
cessamento plásmico de circuitos 
integrados. Os grãos subsequente- 
mente causam efeitos poluidores 
nos circuitos. O mecanismo de for- 
mação de grãos nestes processos 
industriais começa por uma etapa 
de nucleação em que se formam 
“cristalitos”. Depois segue-se uma 
coagulação que dá origem a parti- 
culas maiores. A etapa final de 
crescimento realiza-se à custa de 
deposição na superfície de modo a 
formarem-se os grãos amorftos. 
Como consequência do potencial 
eléctrico na câmara de processa- 
mento os grãos formam nuvens lo- 
calizadas. Dai que o controlo da for- 
mação e localização de grãos em 
plasmas é crucial em muitas tecno- 
logias industriais e requer um co- 
nhecimento não só da física dos 
plasmas básica como também da 
quimica que é relevante em astrofi- 
sica. 


Nos anos 40, Landau demonstrou 
que um sistema de grãos carrega- 
dos, de determinada densidade, 
possui uma energia interna mais 
baixa se este se organizar numa 
estrutura regular em rede. Quando 
a razão entre a energia potencial 
coulombiana e a energia cinética 
dum sistema de particulas carre- 
gadas excede um valor crítico G. 
= 170 forma-se uma rede coulom- 
biana. Dado que o plasma tem um 
efeito de blindagem sobre o 
campo das partículas, a condição 
de solidificação é valida apenas 


quando a distância entre partí- 
culas é menor ou igual à distância 
de Debye. Experiências com sus- 
pensões electrolíticas confirmaram 
esta ideia. Contudo, pensava-se 
que um cristal plásmico só se po- 
deria formar em condições de mi- 
crogravidade, no espaço por 
exemplo. Recentemente, realiza- 
ram-se observações de camadas 
de partículas levitadas no invólu- 
cro de descargas de rádio fre- 
quência de Argon. Por observação 
microscópica e difusão laser pode 
confirmar-se que estas camadas 
encontravam-se organizadas em 
estruturas regulares para dados 
parâmetros da descarga. A título 
de exemplo, observou-se grãos 
com raios de aproximadamente 
10 microns organizados num rede 
coulombiana hexagonal com uma 
distancia entre grãos de cerca de 
880 microns. A carga dos grãos, 
medida pela sua resposta resso- 
nante a um campo periódico ex- 
terno, é da ordem dos -20 x 10º 
cargas electrónicas. Estes siste- 
mas são extremamente interes- 
santes porque permitem o estudo 
de líquidos e sólidos bem como 
transições de fase praticamente a 
olho nu. 

Os desafios imediatos consistem no 
crescimento controlado e caracteri- 
zação dos respectivos diagramas 
de fase e no estudo dos modos pró- 
prios destes sistemas. Estes novos 
cristais diferem dum modo essencial 
dos cristais convencionais por duas 
razões; em primeiro lugar porque as 
cargas dos grãos são variáveis in- 
dependentes e em segundo lugar 
porque estes sistemas são sistemas 
abertos, isto é requerem um influxo 
de energia para manterem a estru- 
tura cristalina. Isto acarreta implica- 
ções tanto práticas com teóricas. 
Do ponto de vista das aplicações, 
estes novos materiais podem apre- 
sentar propriedades uteis até agora 
desconhecidas. Do ponto de vista 
teórico tais cristais exemplificam sis- 
temas simples de estruturas coeren- 
tes, ou seja estruturas do tipo que 
ocorrem nos seres vivos. 


NOVOS MATERIAIS 
PARA APLICAÇÕES 
EM ENGENHARIA 


Rogério Colaço 
Dep” de Engenharia de Materiais, IST 
Original Submetido em: 3 de Fevereiro de 1995 


Por uma questão de sistematiza- 
ção, é comum dividirem-se os ma- 
teriais para aplicações em engenha- 
ria em, pelo menos, três classes: 
polímeros, cerâmicos e metais. 
Destas três classes de materiais, os 
metais foram os que mais condicio- 
naram a evolução das sociedades 
humanas durante os últimos 5000 
anos e, por isso, os antropólogos 
recorrem aos termos Idade do Co- 
bre, Idade do Bronze e Idade do 
Ferro para classificar os estados de 
desenvolvimento das civilizações e 
para fazer a distinção entre a era 
modema e o neolítico (Idade da Pe- 
dra). Desde a Idade do Cobre 
(cerca de 3800 AC) até meados do 
século passado, o domínio das téc- 
nicas da extracção dos metais, da 
fundição e da forjagem, de impor- 
tância estratégica vital, era pertença 
dos mestres artífices que, cuidado- 
samente, transmitiam os conheci- 
mentos acumulados ao longo dos 
anos aos seus melhores aprendi- 
zes. Contudo, os artífices, indepen- 
dentemente da sua perícia, não sa- 
biam porque é que a lâmina da 
espada, quando aquecida até ao ru- 
bro e enfiada numa tina de água fi- 
cava dura, frágil e resistente e, 
quando deixada durante a noite sob 
as brasas da forja, ficava macia e 
facil de deformar, pelo martelo e bi- 
gorna, na manhã seguinte. A revo- 
lução científica do século XVll e a 
revolução industrial do século XVIII 
não se refletiram de imediato sobre 
as técnicas metalúrgicas. Se bem 
que em 1627 Jousse faça as pri- 
meiras observações com carácter 
científico das propriedades dos me- 
tais e, por volta de 1722, Réaumur 
tente relacionar as propriedades de 
um ferro fundido com a estrutura 


que observa com o auxílio de um 
microscópio [1], a verdade é que 
“durante os séculos XVll e XVIII os 
alquimistas não falavam com os fer- 
reiros” como refere G. Liedl [2]. 
E apenas em 1866 que Henry Clif- 
ton Sorby se apercebe que a tex- 
tura da superfície de fractura dos 
aços depende dos tratamentos tér- 
micos a que estes haviam sido sub- 
metidos e da sua composição qui- 
mica exacta. Após polimento e 
contrastação da superfície, Clifton 
Sorby observa pela primeira vez a 
microestrutura metalográfica do aço 
e descreve a perlite. São os traba- 
lhos de Clifton Sorby que estão na 
origem da metalurgia física ou, mais 
genericamente, da ciência dos ma- 
teriais, tal como hoje são entendida- 
das [3] 

Nos cerca de 100 anos que se se- 
guiram a Clifton Sorby até aos nos- 
sos dias, o desenvolvimento de 
uma série de novos instrumentos 
de análise, como os microscópios 
electrónicos de transmissão e de 
varrimento e o difractôómetro de 
raios-X, permitiu aos cientistas 
aperceberem-se da hierarquia de 
niveis estruturais existentes nos 
materiais e correlacionarem essas 
caracteristicas estruturais com a 
maioria das propriedades observa- 
das. O crescente domínio do binó- 
mio estrutura/propriedades e o de- 
senvolvimento das tecnologias de 
produção e de análise fez com que 
o desenvolvimento da ciência e da 
engenharia de materiais se tor- 
nasse um processo de feedback 
positivo: novos materiais permitiram 
a construção de novos equipamen- 
tos, que permitiram “descobrir” no- 
vas características e que “criaram” 
novas necessidades, que conduzi- 
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ram ao aparecimento de novos ma- 
teriais ou de materiais com novas 
caracteristicas. 


À evolução da ciência e da enge- 
nharia de materiais, durante as ulti- 
mas décadas, tem sido de tal forma 
rápida que, hoje em dia, dizer-se 
que que os metais são dúcteis e 
bons condutores, que os cerâmicos 
são frágeis e isolantes ou que os 
polímeros são pouco densos, iso- 
lantes e pouco resistentes começa 
a deixar de ter sentido. Os cerâmi- 
cos supercondutores e semicondu- 
tores, as ligas metálicas leves, os 
metais orgânicos [4] ou os compósi- 
tos (hibridos que reunem as “melho- 
res” características dos materiais 
que os constituem) são alguns dos 
exemplos destes novos materiais 
cujas propriedades não se “encail- 
xam” na classificação tradicional. 

Muitos destes novos materiais tor- 
naram-se indispensáveis ao funcio- 
namento das sociedades modernas» 
e alteraram radicalmente o modo 
de vida das populações. Os cerá- 
micos semicondutores permitiram a 
substituição do tubo de vácuo e ori- 
ginaram o espantoso desenvolvi- 
mento da microelectróônica, as su- 
perligas metálicas, que suportam 
sem perdas significativas das suas 
propriedades mecânicas temperatu- 
ras até 1600ºC, tornaram quase 
banais os voos superssónicos e as 
missões espaciais, os polímeros 
sintéticos mais ou menos sofistica- 
dos que nos rodeiam numa enorme 
diversidade de aplicações e as fi- 
bras ópticas que estão a originar 
uma nova “revolução” da indústria 
de telecomunicações são alguns 
dos (muitos) exemplos possíveis. 
Outros novos materiais, apresen- 
tam propriedades que abrem pers- 
pectivas interessantes prevendo-se 
o alargamento da sua utilização a 
curto/médio prazo. Exemplos disso 
são os vidros metálicos, sintetiza- 
dos pela primeira vez por Pol Du- 
wez no início dos anos sessenta 
[5]. Os vidros metálicos, obtidos 
submetendo algumas ligas metali- 
cas a taxas de arrefecimento entre 
10º e 10'º K/s, apresentam eleva- 


das resistências mecânicas e ca- 
racterísticas magnéticas especiais, 
sendo crescentemente utilizados 
quer como fibras de reforço de 
compósitos quer em aplicações 
que tiram partido do facto do seu 
campo coercivo ser excepcional- 
mente baixo (da ordem de alguns 
milicersteds) [6]. Outro exemplo, 
mais recente e ainda em fase de 
desenvolvimento, são os nanocom- 
pósitos. Estes materiais são obti- 
dos manipulando pós cerâmicos ou 
metálicos com dimensões compre- 
endidas entre 1 e 50 nanômetros e 
apresentam, de uma forma geral, 
uma melhoria significativa das pro- 
priedades mecânicas e magnéticas, 
entre outras, devido à possibilidade 
de “construção” da microestrutura, 
de forma controlada, desde escalas 
quase atômicas [7]. 


Este “Big Bang” de novos materiais 
originou já a proposta de designa- 
ção de Idade dos Novos Materiais 
(Advanced Materials Age) para a 
era em que vivemos [3]. Embora a 
sugestão tenha partido de cientistas 
de materiais, a verdade é que a 
existência deste largo (e crescente) 
espectro de materiais disponíveis, 
introduziu alterações radicais no 
tipo de equipamentos que podem 
ser construidos, no seu design e na 
metodologia de selecção de materi- 
ais para projecto. Uma vez que a 
relação univoca existente entre o 
equipamento a ser fabricado e o 
material utilizado deixou de fazer 
sentido, os engenheiros deparam- 
-Sse, com dois novos problemas: 


* Como seleccionar racionalmente 
o material a utilizar para uma 
nova aplicação”? 

* Quando é que um material tradi- 
cional utilizado numa determi- 
nada aplicação deverá ser subs- 
tiuido por um novo material? 


A resposta a estas duas questões 
torna necessária a introdução do 
conceito de “adequação ao uso” [8]. 
Neste conceito aparecem como fac- 
tores a considerar na selecção de 
materiais não só as suas proprieda- 
des, mas também factores como o 


preço, o processo de fabrico, a du- 
rabilidade, o desempenho em ser- 
viço, a facilidade de reparação, O 
design e um factor de crescente im- 
portância que é o impacto ambiental 
e a facilidade de reprocessamento 
ou reciclagem do material (recorde- 
-se as recentes campanhas publici- 
tárias em torno dos automóveis 
100% recicláveis ou dos detergen- 
tes “verdes”). À introdução deste 
conceito implica que a selecção de 
materiais se torne necessariamente 
uma tarefa multidisciplinar, de modo 
a que os diversos aspectos da es- 
colha possam ser levados em con- 
sideração. 


Esta necessidade de perspectivar 
um determinado material conside- 
rando as suas múltiplas dimensões 
(propriedades, processo de fabrico, 
preço, etc.) tem como consequên- 
cia que alguns dos novos materiais 
sejam “colocados na prateleira” 
após o entusiasmo que se segue à 
sua descoberta. Dois exemplos 
disso são as ligas com memória de 
forma e os cerâmicos supercondu- 
tores de alta temperatura crítica 
(CSC). O efeito de memória de 
forma foi pela primeira vez obser- 
vado numa liga TiNi em meados 
dos anos 50 [9], tendo sido expli- 
cado cerca de 10 anos mais tarde, 
por Lange e Zijderveld [10], com 
base na existência de transforma- 
ções de fase reversíveis e reprodu- 
tiveis (transformações termoeslásti- 
cas) que permitiam à liga recuperar 
a sua forma inicial após ter sido de- 
formadas plasticamente, mediante 
um simples aquecimento. Nos anos 
que se seguiram ao trabalho de 
Lange e Zijderveld surgiram inúme- 
ros resultados de investigação so- 
bre este tipo de ligas, tendo sido 
sugeridas inúmeras possibilidades 
para a sua utilização. Hoje em dia 
a utilização das ligas com memória 
de forma limita-se a pouco mais do 
que o fabrico de pequenas curiosi- 
dades como “a-colher-que-se-endi- 
reita-quando-introduzida-na-chá- 
vena-de-chá”. A razão para tal 
prende-se essencialmente com o 
facto das vantagens comparativas 
das ligas com memória de forma 


ser pouco significativa face a mate- 
riais tradicionais, na maioria das 
possiveis aplicações, e o seu custo 
ser relativamente elevado. O caso 
dos cerâmicos supercondutores é 
um pouco diferente. A descoberta 
da supercondutividade para tempe- 
raturas da ordem dos a 35 K numa 
liga cerâmica Ba-La-Cu-O por Bed- 
norz e Muller em 1986 [11], foi tal- 
vez das descobertas que maior agi- 
tação provocou na comunidade 
científica nos últimos anos, tendo 
valido o prémio Nobel a estes dois 
investigadores. Uma vez que se 
previa que os supercondutores re- 
volucionassem as tecnologias da 
electrônica, dos transportes (devido 
ao efeito de exclusão de campo 
magnético aplicado que os CSC 
apresentam é possível construir 
comboios que “levitam sobre os 
carris”), dos sistemas de potência, 
etc., foram feitos avultados investi- 
mentos em I&D nestes materiais e, 
em 1987, a temperatura crítica dos 
CSC já tinha ultrapassado os 90 K 
(numa liga Y-Ba-Cu-O que apre- 
sentava uma temperatura crítica 
93 K). Hoje, cerca de 10 após a 
descoberta de Bednorz e Muller, o 
entusiasmo em volta dos CSC “ar- 
refeceu” um pouco. As dificuldades 
de produção dos CSC e a sua ex- 
trema fragilidade, já para não men- 
cionar o facto da temperatura crí- 
tica ainda não ter ultrapassado os 
125 K, colocam importantes obstá- 
culos à sua utilização generalizada 
a curto prazo. 

Por outro lado é preciso recordar 
que os desenvolvimentos da ciência 
e engenharia de materiais não tem 
como produto apenas o apareci- 
mento de novos materiais: tem per- 
mitido também alterações das pro- 
priedades dos materiais tradicionais 
mediante o desenvolvimento de no- 
vas tecnologias. Por exemplo, o de- 
senvolvimento das tecnologias de 
superfícies (tratamentos superficiais 
por laser, CVD, PVD, nitroração ió- 
nica, etc.) permitiu, entre outras coi- 
sas, aumentos substanciais da re- 
sistências ao desgaste, à corrosão 
e à oxidação de materiais tradicio- 
nais, sem que o processo ou o 
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custo de fabrico sejam significativa- 
mente alterados. Este factor pode 
ser decisivo e fazer com que um 
novo material seja preterido por um 
material tradicional numa determi- 
nada aplicação: Uma vez que as 
tecnologias dos materiais tradicio- 
nais são, de um modo geral, melhor 
dominadas, melhorar as proprieda- 
des de um material já utilizado im- 
plica, normalmente, um menor es- 
forço de reconversão do processo 
de fabrico, do que a sua substitui- 
ção por um novo material. 


Em conclusão poderá dizer-se que 
o desenvolvimento da ciência e da 
engenharia de materiais fez com 
que a oferta de novos materiais e 
de materiais com novas proprieda- 
des aumentasse extraordinaria- 
mente. Sendo expectável que este 
“Big Bang' de materiais continue, 
torna-se necessário que algum es- 
forço de I&D seja canalizado de 
modo a assimilar esta tão grande 
diversidade, quer em termos de se- 
lecção racional de materiais quer 
em termos de desenvolvimento de 
tecnologias de produção de novos 
materiais economica e ambiental- 
mente viáveis. 
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RESUMO 


A presente série de dois arti- 
gos aborda o problema da in- 
trusão e deposição de sal ma- 
rinho nas áreas terrestres, um 
tópico de contaminação ambi- 
ental que aparece frequente- 
mente desvalorizado face a 
outras cargas naturais ou an- 
tropogénicas. Nesta primeira 
parte, faz-se uma breve revi- 
são das consequências tec- 
nico-econômicas da salinidade 
atmosférica, centrada no res- 
pectivo impacto sobre os ma- 
teriais de Engenharia e sobre 
a generalidade das infra-estru- 
turas produtivas; na segunda 
parte, descrevem-se as primei- 
ras etapas do processo de se- 
lecção e avaliação de epifitos 
criptogâmicos para eventual 
utilização como biomonitores 
da densidade de deposição 
salina na litosfera. 


ABSTRACT 


Deleterious effects of airborne 
salinity are currently observed 
on man-made structures, vege- 
table organisms (crops and na- 
tives) and soil/water resources; 
to state only a few, such ef- 
fects include accelerated de- 
gradation of materials and vari- 
ous contributions to land 
desertification, namely through 
salt-assisted deforestation, 
aquifer contamination and di- 
rect soil salinizaton. This article 
outlines some actual and po- 
tential hazards associated with 
salt-laden atmospheres and 
briefly reviews the economic 
impact of chloride-induced at- 
mospheric corrosion. 


Palavras chave: biomonitores, contaminação natural, corrosividade at- 


mosférica, impacto salino, materiais, salinidade tropos- 


férica. 
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RISCO SALINO 


Com um fluxo troposférico estimado 
em 10º(-10'º) t/'ano [1,2], o sal mari- 
nho constitui-se no maior conjunto 
de material fragmentário em circula- 
ção na atmosfera, atingindo quase 
40% do total admitido para todas as 
partículas — primárias e secundárias 
— de todas as origens — naturais e 
artificiais — cuja emissão devera ul- 
trapassar 2600 Mt/ano [3,4]; parte 
significativa do sal acaba por ser 
dispersada sobre as massas terres- 
tres mediante processos de deposi- 
ção seca e húmida, sendo que os 
primeiros poderão representar cerca 
de 3/4 do fluxo reentrante total [5]. 
Indispensável ao equilíbrio ambien- 
tal dos ecossistemas costeiros e 
vector de biodiversidade e zonagem 
nos mesmos [6-14], a salinidade at- 
mosférica não deixa de ser um fac- 
tor de risco para os materiais (so- 
bretudo metálicos), a fito-sanidade 
(vegetação e culturas) e a qualidade 
dos solos e/ou recursos hídricos. 


As orlas costeiras têm um peso de- 
cisivo na geografia humana e eco- 
nómica a nível mundial: aí se 
encontram fixados importantes seg- 
mentos demográficos, e aí se en- 
contra implantada grande parte do 
aparelho produtivo e do equipa- 
mento social; acresce que, embora 
a respectiva agressividade seja va- 
riável, as atmosferas marinhas são, 
em regra, bastante severas para 


com as estruturas industriais-urba- 
nas. Estes factos deveriam levar, 
por si sós, a que o impacto salino 
viesse a ser sistemática e conveni- 
entemente considerado nos projec- 
tos envolvendo tais áreas, mas não 
só: a influência maritima não se 
circunscreve às regiões periféricas 
e pode estender-se muito para o 
interior [14], chegando autores de 
vários quadrantes a referir cente- 
nas [15,16] ou mesmo milhares de 
quilómetros [17], isto para não falar 
da especificidade dos sistemas in- 
sulares [18-20]. Nestas condições, 
a avaliação e previsão do risco sa- 
lino afigura-se determinante para 
uma adequada selecção de materi- 
ais e plantas, em particular, e para 
qualquer estratégia integrada de 
protecção, reabilitação e harmoni- 
zação de áreas sensíveis, donde, 
no geral, para uma mais equili- 
brada gestão de espaços/recursos 
naturais e um melhor ordenamento 
do território. 


IMPACTO SALINO 


Uma pálida imagem da relevância 
(negativa) da deposição salina na li- 
tosfera poderá ser obtida mediante 
consideração, por exemplo, do res- 
pectivo impacto sobre os materiais 
metálicos º, uma vez que a degra- 
dação atmosférica responde por 
mais de metade dos custos totais 
da Corrosão [21]. Os números dis- 


3 Muitos outros exemplos poderiam ser considerados: ainda recentemente, através 
do respectivo programa específico de investigação e desenvolvimento ambiental 
(ENVIRONMENT, 1991-94), a União Europeia indiciava a salinidade como um dos 
factores abióticos de desertificação na área mediterrânica. 
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poniveis são, no mínimo, impressio- 
nantes: desde a elaboração do es- 
tudo NBS-Battelle [22] que nos pai- 
ses desenvolvidos se admite, em 
média, 3.5% do respectivo PNB 
para valor total daqueles custos... 
Saliente-se que, a despeito de algu- 
mas estimativas mais conservado- 
ras, tanto os EUA [23] como a CEE 
[24] ainda aceitam o estimador glo- 
bal (4.2% do PNB) decorrente do re- 
ferido estudo. Em Portugal, a situa- 
ção não se encontra avaliada, muito 
embora não deva ser mais favorá- 
vel, isto porque: 


|) quer no Continente (extensão 
da costa, regime de ventos) 
quer, obviamente, nas Ilhas, 
faz-se sentir uma forte influên- 
cia oceânica; 

li) a assimetria no desenvolvi- 
mento litoral-interior leva a que 
a geografia socio-econômica 
seja predominantemente cos- 
teira; 

ii) nem a Administração Pública 
nem as Empresas estarão por- 
ventura sensibilizadas para a 
potencial redução nas despesas 
públicas/privadas resultante da 
implementação criteriosa de tec- 
nologia anticorrosiva, redução 
essa que se crê poder atingir 
40% da factura total de Corro- 
são [25]. 


Embora as perdas atmosféricas não 
se encontrem desagregadas por 
factores de risco, parece evidente 
que a parcela imputável à indução 
salina devera responder por parte 
significativa da factura: não é por 
acaso que a taxa de deposição de 
cloreto é uma das três variaveis pri- 
márias de avaliação e classificação 
da corrosividade atmosférica — veja- 
se [26], por exemplo — com as agra- 
vantes de o sal possuir: 


) uma fonte deslocalizada, per- 
manente e de potência incontro- 
lável; 


ii) uma quase auto-suficiência na 
criação/manutenção de filmes 
de electrólito, decorrente das 
características higroscópicas e 
deliquescentes do aerosol mari- 
nho [27-29], o que remete para 


um abaixamento considerável 
do limiar de risco (humidade crí- 
tica) para corrosão significativa 
[30-32]; 

ii) uma agressividade especifica 
para certos materiais e/ou estru- 
turas, cuja deterioração coloca 
sérios problemas tecnico-econó- 
micos ao nível qualitativo e/ou 
quantitativo: veja-se, por exem- 
plo, o caso dos aços patináveis 
— metais muito interessantes do 
ponto de vista da Engenharia/ 
[Economia/Ecologia (alta resis- 
tência, boa formabilidade, baixa 
liga, auto-protecção, reduzida 
manutenção), mas com fraco 
desempenho em atmosferas 
cloretadas [33-35] — ou do be- 
tão armado — provavelmente o 
material de Engenharia com 
custos mais elevados em ter- 
mos de Corrosão [25]; 


iv) uma série de preocupantes si- 
nergias com poluentes de enxo- 
fre e azoto, e com a precipita- 
ção ácida visivel ou oculta 
[36-50]. 


AVALIAÇÃO DO FLUXO SALINO 


Apesar dos seus efeitos directos e 
indirectos, a contaminação salina 
tem merecido um tratamento inferior 
ao que vem sendo reservado a ou- 
tras formas de poluição (gases, po- 
eiras): esta situação — que não pa- 
rece aceitável ainda que se invoque 
tratar-se de um vector natural — po- 
derá ter origem na fraca elastici- 
dade dos modelos empíricos e/ou 
deterministicos de dispersão, e em 
limitações tecnico-econômicas das 
redes clássicas de amostragem. 
Acresce que, no primeiro caso, es- 
taríiamos mesmo perante um falso 
argumento: risco natural não signi- 
fica, de todo, risco menor, a não ser 
que se considerem igualmente me- 
nores os riscos sísmico e vulcânico, 
o potencial de auto-ignição dos 
ecossistemas, ou a emissão não- 
antropogênica de compostos de en- 
xofre e azoto [51-56]: já o(s) se- 
gundo(s) aspecto(s) se prende(m) 
mais com insuficiências materiais 


na modelização da densidade de 
deposição salina, agravadas pela 
extrema variabilidade inerente à 
produção, transporte e distribuição 
do aerosol marinho. 


O emprego de plantas como indica- 
dores de salinidade poderá levar a 
que, no futuro, aquelas redes ve- 
nham a ser operadas de um modo 
flexivel, descentralizado e ecológico. 
Com efeito, a substituição dos cap- 
tadores físicos por biomonitores 
perspectiva apreciáveis vantagens, 
desde a possibilidade de refina- 
mento continuo da malha de amos- 
tragem e de replicação amostral 
praticamente ilimitada, até à dimi- 
nuição nos encargos de instalação, 
manutenção e operação relativa- 
mente às estações tradicionais, elas 
próprias algo intrusivas em termos 
paisagísticos. De resto, as óbvias li- 
mitações no mapeamento salino 
com base nestas últimas poderão 
estar na origem do fraco desempe- 
nho que vem sendo atribuído aos 
modelos disponíveis para a intru- 
são/dispersão de sal, os quais se 
ressentem da escassez de resulta- 
dos de campo mesmo nos casos 
em que a logistica não parece ser 
excessivamente pesada [57]. 


Na segunda parte deste trabalho, 
serão apresentadas e discutidas as 
primeiras fases do estudo que vem 
sendo efectuado sobre epífitos crip- 
togâmicos, enquanto candidatos a 
indicadores de intrusão salina at- 
mosférica e monitores da respectiva 
densidade de deposição terrestre. 
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ABSTRACT 


Following the miniaturization of 
electronic devices, thin film so- 
lid state batteries are receiving 
some new attention, motivated 
by the desire for a battery that 
can be fully integrated with mi- 
crocircuits. Lithium microbatte- 
res have been the most stu- 
died; while some of them are 
the replica of “classical” batte- 
ries using a liquid electrolyte, 
e.g the Li/Tis. system, others 
use new materials. Very inte- 
resting results have been ob- 
taineo. 


Because of the very high con- 
ductivity of silver conducting 
electrolytes silver microbatte- 
ries have also been extensi- 
vely studied and significant re- 
sults obtained with a primary 
cell design. 


Other systems using solid 
state ionic conductors have 
also been considered; they are 
characterized by a low availa- 
ble voltage. 


THIN FILM SOLID 


STATE MICROBATTERIES 


L.P.5. Araújo, M.J.C. Plancha and C.A.C. Sequeira 
Original Submetido em: 12 de Outubro de 1994 
Instituto Superior Técnico, Lisboa 


1. INTRODUCTION 


Thermodynamically, an electrochemical e.m.f. system (a so-called battery) is 
generated if an electrolyte is sandwiched between two electrode materials 
with different chemical potentials. Further, if a constant supply of ions can be 
maintained and transported through the electrolyte, it will deliver current 
when connected across a load resistance. 


Two different kinds of batteries are used, primary (non-rechargeable) and se- 
condary (rechargeable); they comprise liquid or solid electrodes and elec- 
trolytes. 


Solid state batteries have been studied for about 30 years. Their advantages 
compared to those using a liquid electrolyte are now well known: ease of uti- 
lization, resistance to shocks and vibrations, absence of possible pollution 
due to liquid electrolyte, thermal stability, absence of self-discharge (owing to 
the very low electronic conductivity of the solid electrolytes generally used) 
and possible miniaturization [1]. 

The current densities obtained with these batteries were generally low be- 
cause of the high internal resistance of the cell, which mainly results from 
the low ionic conductivity of the electrolyte [2, 3]. 

Up to now, only three types of solid state button cells have been developed. 
These are lithium batteries in which the electrolyte is Lil [4], a composite Lil- 
Al,O, [5] and a sulphur glass [6]. One method of reducing the internal resis- 
tance is to use an electrolyte in the form of a thin film. 


This was the first reason for the use of thin film techniques in battery techno- 
logy. Presently and following the miniaturization of electronic devices, thin film 
solid state batteries are receiving some attention, motivated in part by the de- 
sire for a battery that can be fully integrated with microcircuits. This concept 
involves a power source and circuit constructed on the same chip or subs- 
trate, employing either a single integrated battery or a distributed power 
source. The operational acceptability of such integrated devices is contingent 
on the low power requirements of the new circuits and on the high energy 
density of the emerging lithium batteries. However, thin film batteries have low 
capacities. Rechargeability is an essential condition for technical feasibility. 
Another advantage of thin film cells is the ability to make large area batteries 
easily so that low current densities are acceptable. 


In addition, thin film technologies can allow for the preparation of new mate- 
rials. Thin film technologies are divided into two main categories: the first 
concerns physical deposition methods, including vacuum deposition, sputte- 
ring and ion plating; the second concerns chemical deposition methods, 
mainly chemical vapour deposition (CVD). It is well known that, using physi- 
cal deposition methods such as evaporation, the “vapour” composition can 
differ strongly from that of the starting material [7]. The condensation of 
these vapours can lead to new materials, perhaps providing better perfor- 
mances. Consequently it is absolutely necessary to characterize the thin 
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films obtained, bearing in mind that two main difficulties remain: the film 
thickness and the very small amount of total material within the thin films. 
Table 1 gives a view of the main techniques used to characterize thin films. 
This paper will try to give an overview of thin film microbatteries. The term 
microbatteries is used to describe devices constructed using thin film techno- 
logies, and typically each part of the cell — both the electrodes and the elec- 
trolyte — comprise thin films of solids a few micrometres thick. 


Table 1 

Main techniques used to characterize thin films (from ref. 8) 

Instrumental method Acronymi(s) Morphology Chemical Crystallographic Electrical 

determination analysis structure and mapping 

mechanical 
properties 

Auger electron AES X 

spectroscopy 

Electron-beam-induced EBIC X 

current microscopy 

Laser reflectance LR X 

Neutron activation analysis NAA X 

Normarski interference — X 

contrast optical microscopy 

Rutherford backscatteling RBS X (X) 

spectroscopy 

Scanning electron SEM X (X) (X) 

microscopy 

Secondary ion mass SIMS X 

spectroscopy 

Transmission electron TED X 

diffraction 

Transmission electron TEM X (X) X (X) 

microscopy 

Voltage contrast VC X 

microscopy 

X-ray diffraction XRD X 

X-ray emission XES X 

spectroscopy 

X-ray fluorescence XRF X 

X-ray photoelectron XPS, ESCA X 

spectroscopy 


2. LITHIUM MICROBATTERIES 


Lithium is a very attractive negative electrode material because of its very 


high mass capacity (3885 A h Kg'), but its reactivity is very high and res- 


tricts the choice of electrolyte. 

One of the first lithium primary microbatteries was proposed by Liang et al. 
[9, 10]. A Lil Lil Agi battery was formed by the following procedure. Agl was 
sprayed on a silver foil and run through a furnace to melt the Agl film on the 
silver foil. Lil electrolyte was vacuum deposited. To prevent internal shorting, 
the thickness of pinhole-free Lil film must be greater than 5um. The lithium 
was then evaporated on the electrolyte. The open-circuit voltage (OCV) of 
the cell of 2.1 V agreed with the value of 2.09V calculated from the free- 
energy change for the reaction Li + Agl => Lil + Ag, and it confirmed this to 
be the reaction of the cell. The formation of Lil as a reaction product implies 
that the impedance of the cell increased sharply during discharge [9, 10]. 
The polarization curve of such a battery is given in Fig. 1. 
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Fig. 1: Typical polarization curve of Li [Lil | Agl cell (from ref. 9) 


An improvement of this system was proposed by Feuillade [11]. One of the 
failures of Liang's battery was the thickness of the Lil thin film of electrolyte. 
As previously shown, to suppress the defects (pinholes, clefts) of the evapo- 
rated thin film it was necessary to increase its thickness to 5 um, thus incre- 
asing the internal resistance of the cell. 


Feuillade proposed chemical oxidation of the lithium by a solution of iodine in 
thionyl chloride to produce a 1 um thin film of Lil. Another method proposed 
by the author is electrochemical oxidation of the lithium. The iodization solu- 
tion was a solution of ammonium tetrabutyl iodide in thionyl chloride. Lithium 
was the anode and platinum a counterelectrode. Initially galvanostatic (until a 
thickness of Lil of about 1 um), the reaction was then potentiostatic to incre- 
ase the compacity of the film without variation of the thickness [1 1). 


Another way to prepare the Lil thin film is to add a small amount of a binder 
to a propylene carbonate solution of Lil. The paste is deposited on Ag/Agl 
thin films. The solvent is evaporated, a foil of lithium is added and the entire 
system is then pressed [12 - 14]. 

In the same system, Rea proposed the substitution of bismuth for Agl be- 
cause of its high electronic and lithium ion conductivities as well as its relati- 
vely high voltage potential with a lithium anode [15]. 

The Li|Lil|Bi cell was entirely constructed by sequential vacuum evapora- 
tion. Such a deposition method provided for uniformity in layer thickness and 
increased ease of fabrication. The main advantage of such a cell is the use 
of bismuth because of the large interspatial openings in its crystalline struc- 
ture capable of easily accommodating anode metal cations; furthermore, bis- 
muth may function as its own current collector. The use of arsenic, antimony, 
lgad, tin, selenium and tellurium instead of bismuth is also possible. It seems 
that there is no change in the crystal structure of the cathode materials, gi- 
ving the further advantage of providing rechargeable cell systems. 


The OCYV of the cell is 0.83 V and a typical current density of 10-100 uA em? 
can be easily obtained, but the reversibility seems to be limited to about 10 
cycles [15]. 

In 1983 Levasseur et al. proposed a lithium battery with a positive electrode 
made from a thin pressed pellet of TiS. [16]. An evaporated borate-based li- 
thium glass was used as electrolyte and an evaporated thin film of lithium 
served as the negative electrode. Such a type of battery could deliver 10 uA 
for 30h with an average voltage of 2 V, but the system had poor rechargea- 


37 


38 


bility [16]. The main problem was the orientation of the 2D shape grains of 
TiS., perpendicularly to the lithium ion motion in the cell (crystallo-graphic c 
axis of TiS, perpendicular to the substrate). 

During the same period Kanehori et al. proposed a solution to this difficulty 
[17-19]. A TiS. film was fabricated by CVD at a substrate temperature of 450 
“C using a gas composed of TiCl, and H,S diluted with helium. The micros- 
tructure of the 1-3 um thick film of TiS., consists of small narrow plate-like 
crystals having their crystallographic c axis parallel to the substrate plane. 
An amorphous lithium phosphosilicate electrolyte and a lithium negative elec- 
trode were evaporated on this layer [20]. This arrangement allowed a good 
reversible intercalation of the lithium in TiS, and can give a capacity of more 
than 100 uA h cm*º. Figure 2 shows discharge and charge curves of this 
type of thin film battery. Depending on the current density and the depth of 
discharge, 100-1500 cycles were possible [20]. 
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Fig. 2: Discharge and charge curves of thin film cells (from ref. 17): O discharge; A charge; 
cell A, TiS, thickness tum; cell B, TiS, thickness 3.7um. 


Another method to obtain TiS, is by codeposition of nickel and TiS, grains on 
a nickel substrate. The electrochemical reduction of dispersed TiS., powder in 
a non-aqueous solution of nickel chloride produced a 10-15 um thick film in 
which the TiS, grains had any direction and were embedded in nickel [21]. 
More recently Jourdaine et al. proposed a microbattery using a mixed con- 
ducting glass obtained in the system V,O.-P,O, or V,O.-TeO, as the positive 
electrode [22]. Successive thermal evaporation depositions of this electrode, 
of lithium borophosphate glass as electrolyte and of lithium as the negative 
electrode formed a microbattery able to deliver a few uA cm? at an average 
voltage of 1.3 V [22]. 

Meunier et al. used r.f sputtering to obtain a new titanium oxysulphide thin 
film material on which a glass (B,O0,-Li,O-Li,SO,) was then deposited with 
the same technique. A thermal evaporation of lithium made the microbattery 
complete [23]. A current density of up to 60 uA cm“ at an average voltage of 


2 V could be easily obtained. This system is highly reversible and more than 
100 cycles have been obtained without trouble [23, 24]. 


Relatively recently many papers have reported the use of thin films of Che- 
vrel phases Cu,Mo,Sa. (2<y<4, 7.8<8-z<8) as positive electrodes. These thin 
films are prepared by the technique of painting the powder on a Cu,S layer 
and subsequent reaction with H.S-H. gas at 530ºC [25-27]. Chevrel phases 
can also be prepared by CVD [28]. Up to 18 lithium ions per unit formula can 
be incorporated, giving a theoretical energy density of 740 W h Kg" at 0.5 
mA cm? current density [26]. However, only two of these are reversible. 
There is no doubt that these films will be used in microbatteries. 


3. SILVER MICROBATTERIES 


Silver halides and rubidium silver iodide provide a very high Ag* conductivity. 
RbAg.l; exhibits a conductivity of 27 Q' em" at 25ºC, which is the highest 
value for all solid electrolytes at room temperature. Although the use of such 
electrolytes in high power solid state batteries is not realistic because of the 
price of silver, their utilization as thin films in microbatteries is a possibility. 
Thin films of silver compounds as the electrolyte have been reviewed by 
Kennedy [29]. 

Vouros and Masters first reported a thin film silver battery [30, 31]. A strip of 
sputtered platinum was followed by the deposition of a film of silver (0.1 um). 
The system was then placed in an iodine chamber to ensure the conversion 
of most of the silver to Agl. A second layer of Agl was deposited by evapora- 
tion over the ultrathin iodized silver film which acted as the nucleus for the 
formation of a uniform and void-free deposit of Agl. The silver electrode was 
then deposited by evaporation. The measured OCYV (0.55 V) is consistent 
with a Pt|Agl|Ag concentration cell. A very low current density (less than 
1 uA em?) was available mainly because of the very high internal resistance 
of the cell. 

Many other Ag' conductors sandwiched between silver and another metal 
are reported in the literature. According to the activity of the silver in these 
other metals, the voltage of the cells is low (less than 1 V) and the current 
density never exceeds a few uA em* [32-35]. 


Takahashi and Yamamoto have patented an Ag | RbAgsl; | Te cell giving inte- 
resting results [36]. Six of these cells were stacked entirely by vacuum eva- 
poration. At 25º€ the OCV was 1.2 V. At 5 mA the operating voltage was 
0.9 V [36]. 


Minami reported other silver conductors used in batteries and electrochromic 
devices but without better performances [37]. 


Another kind of silver battery was proposed by Chandra et al. [38]. The elec- 
trochemical chain Ag |NH,Agals | (C+RbAg,l.+KI,) was used. The electrolyte 
was obtained directly on the silver anode by electrolysing an aqueous solu- 
tion of NH,l. The cathode used was a pelletized KI. +C+RbAg,l; mixture. The 
chelate KI(Kl+l,<>Kl,) was used as a source of iodine and was obtained by 
heating the stoichiometric mixture at 120ºC for 36 h. 

The OCV obtained (0.6 V) is lower than that obtained with the classical 
Ag | Ag: ||, configuration (0.69 V). The electrolyte in this case is not a pure io- 
nic conductor. In spite of this disadvantage, several hundred uA em? could 
be obtained at 0.5 V. 


4. LEAD MICROBATTERIES 


Kennedy and Hunter were the first to report a cell using lead metal as the 
negative electrode, PbF, as the electrolyte and CuF, as the positive elec- 
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trode [39]. Because of the high resistivity of CuF., it was codeposited with 
PbF. to form a more conductive mixed layer. The reactions 


[=| Pb+2F”" — PhF,;+2e 


[+] CuF,+2e > Cu+2F 


lead to a voltage of 0.7 V. A current density of up to 10 uA cm? could be ob- 
tained. This cell was not rechargeable. Other works reported microbatteries 
using PbF, (or derived structures) or other F' conductors as the electrolyte 
[40, 41]. Lead was always the negative electrode and BiF, often the positive 
electrode. In these cases the system was rechargeable. Current densities of 
40 uA cm“ at 0.2 V were obtained [40]. 


The use of an evaporated mixture of PbF. and LiF as the electrolyte was 
also mentioned [42]. 


5. COPPER MICROBATTERIES 


Although copper ion conducting materials are well known, applications in thin 
film microbatteries are very rare. 


One example is a battery with bimetal electrodes sandwiching a Cu' conduce- 
tor. The latter was prepared by an electrolytic deposition method giving a Cul- 
-CuCI mixture on a copper anode surface. Many other metals could be used 
on the other side. In all cases a low voltage was obtained, between 0.26 V 
(chromium) and 1.3 V (magnesium) [43]. In spite of the very low conductivity 
of the electrolyte (10º O! em! at 25ºC), a few uA cm were obtained. 


6. THIN FILM BATTERIES USING DIVALENT POLYMER ELECTROLYTES 


lonically conducting polymers consisting of a complex between a metal salt 
and poly (ethylene oxide) (PEO) offer a number of advantages over other 
solid electrolytes, with regard to their possible use in advanced solid-state 
electrochemical devices. Interesting results on these electrolytes have been 


reported by several workers [44-46]. 


In our Laboratory we have prepared and studied polymeric complexes of 
poly(ethylene oxide) (PEO) with salts of zinc, copper and silver [47]. It was 
concluded that some of the studied compounds are good solid electrolyte ma- 
terials in view of their high ionic conductivity, their low electronic conductivity, 
and their chemical compatibility with high energy density electrode materials. 


Experiments have been carried out on cells utilising electrolytes, the conducti- 
vities of which are at the lower limit of values acceptable for batteries for pre- 
sent day use. The choice of type and amount of salt used in these electrolyte 
films was dictated by a systematic study of factors affecting conductivity, the 
results of which have been published elsewhere [3, 48]. The performance of 
the cells was further limited by the use of a simple casting technique to form 
the composite (polymeric electrolyte matrix/active constituent) cathode; it has 
since been shown that significant improvements result ffom modifying the 
casting technique [49, 50]. Despite these limitations, the currents drawn were 
modest but adequate for certain practical devices. An optimised choice of 
anion and stoichiometry for a given divalent cation species can lead to an im- 
provement of two orders of magnitude in conductivity [51], and a concomitant 
improvement in cell performance would be anticipated. 

Some examples of preliminary cells that have been prepared and subjected 
to discharge-charge cycles and other dc measurements are as follows. 


a) Zn(PEO),nClyPolyaniline test cells: The test cells exhibited an open-cir- 
cuit voltage of 1.3 V and a short-circuit current of 80uA (= 80 uA cm?) at 


25 + 1ºC. The test cells were discharged and charged, between 1.3 V and 
0.5 V, at a constant discharge load of 100 Kº and a constant charge cur- 
rent of 10 LA, respectively, for more than 300 hours. The cycles were smo- 
oth and polarization was noted only after the tenth discharge-charge cycle. 


£n(PEO) ,nCl//Methyipolyaniline test cells: The test cells exhibited an 
open-circuit voltage of |.IV and a short-circuit current of 70 uA. The first 
discharge curve was smooth as in the case of the Zn/polyaniline cell, but 
the recharge curve was accompanied by gross polarization and the cell 
was not able to sustain its charge. 


£n(PEO)2nCIsSPEO, TiS, test cells: The Zn/50% TiS., 40% PEO, 5%C 
test cells exhibited an open-circuit voltage of 0.5V and a short-circuit cur- 
rent of 45 uA. Severe polarization losses were noted during the first dis- 
charge and we could not cycle the cells with a minimum of efficiency. 
The cells exhibited a very poor charge retention. 

Zn(PEO)AnCl/ MnO, test cells: The test cells exhibited an open-circuit 
voltage of 1.4V and a short-circuit current of 90 uA, similar to those of the 
Zn/polyaniline cells. Under similar discharge-charge conditions (100 KO, 
10uA, 25+1ºC) the polarization loss of the Zn/MnO, test cells was consi- 
derably less than that of the Zn/polyaniline cells. Figure 3 shows a typical 
discharge curve for one of these test cells at 20ºC. 


CFLL AS 50 "E 


LOAD « 105 E ohm 


Fig. 3: Discharge curve of a laboratory cell Zn | (PEO)ZnCl; | MnO.. 


e) Zn/(PEO),NH,CI.22nCl,/MnO, test cells: The test cells exhibited an open- 


circuit voltage of 1.5V and a short-circuit current of 1.5 mA. The voltage 
drop accompanying the discharges under constant load of 100 KO was 
insignificant, and the rechargeability of the system was demonstrated by 
more than ten discharge-charge cycles performed without serious polari- 
zation effects. 

Zn/(PEO),.KI.22nCl/MnO, test cells: The test cells exhibited an open-cir- 
cuit voltage of 1.1V and a short-circuit current of 1.0 mA. Almost comple- 
tely discharged cells (cell voltage from 1.1V to 0.30V) could not be re- 
charged, and discharge-charge cycles performed between 1.1V and 0.7V 
showed significant polarization, even for very slow cycles. 


An examination of ceils a) to f) revealed that only three systems were reaso- 
nably well rechargeable (Table 2). Furthermore both the Zn/MnO, systems 
operate at higher efficiencies than the Zn/polyaniline system. From the 
energy availabilty the (PEO),ZnCl, is less desirable than the 
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(PEO),.NH,CI.22nCl,. Accordingly, among cells a) to ), we concluded that 
the Zn/(PEO),.NH,CI.2ZnCl,/MnO, is the most suitable system for the deve- 
lopment of practical storage batteries. 


Table 2 
Thin film rechargeable cells 


Zn | (PEO),. ZnCl, | Polyaniline 
Zn | (PEO), ZnCL, | MnO, 
Zn | (PEO), NH,CI.2ZnCI, | MnO, 


Cell cycling consisted of constant load discharge (= 100 kQ) / constant current charge (= 
10 uA) cycles: 

the cells could be cycled 1-10 times, but loss of 

efficiency was observed even for the 1 st cycle. 


After the preliminary preparation and testing of cells a) to f), more adequate 
attention was given to the formation of the composite cathode. The practical 
devices obtained gave the following results. 


9) Mg'PEO,.: Mg(SCN).(45% TiS., 50% PEO,s: Mg(SCN),, 5% acetylene 
black) test cells: The test cells exhibited an open circuit voltage of 1.7 V at 
30 + 1ºC. The test cells were discharged at a constant load at 1 ml until 
apparent exhaustion of active cathode material had occurred. À current of 
1 uA could be maintained for 1000 hours and a slightly lower current 
could be sustained for a further 2000 hours or so. Additionally, test cells 
were used to supply the power to an LCD clock (requiring >5 uA at >1.2 
V) for periods of two weeks at ambient temperatures as low as 15º€. 


Table 3 compares the open circuit voltages and short circuit currents of the 
Mg cell with the data for the Zn cells. There is a significant improvement of 
the Mg cell. 


h) MPEO,s:MX/(45% TiS, or V5O,, 50%PEO,s:MX., 5% acetylene black) 
test cells: These test cells included Mg, Mg(SCN).: Mg(CIO,),; Zn, 
ên(CIO,), anode lelectrolyte combinations. 


nt ed voltages and short circuit currents of some thin film batteries 
Cell Voc, volt Isc, LÃ 

Zn | (PEO),.ZnCI. | Polyaniline RE Bo 

Zn | (PEO),»ZnCI, | Methylpolyaniline 141 70 

Zn | (PEO);ZnCI, | PEO, TiS. 45 

Zn | (PEO)aZnCI, | MnO; 1.4 90 

Zn | (PEO),NH,CI.2ZnCI, | MnO. 1.5 1,500 

Zn | (PEO),KI.2ZnCI, | MnO, 1.1 1,000 

Mg | (PEO),.-Mg(SCN), | PEO, TIS, 1.7 2,500 


For Mg/Mg(SCN),'TiS., a current of 0.5 LA could be sustained for 400 hours. 
For Mg/Mg(CIO,)./V;O,« cells, a current of 1 uA could be maintained for the 


same time. For Zn/Zn(CIO,)»/TiS., a current of 0.5 uA could be sustained for 
300 hours, but if TiS, were replaced by V,O,, the same current could only be 
maintained for 100 hours. 


) M/PEO,:M(SCN)/(45% TiS, or V5O,;, 50% PEO,s: MX., 5% acetylene 
black) test cell: Here the anode metal differed from the divalent cation in 
the electrolyte. For test cells involving a Zn anode with a PEOs: Mg(SCN), 
electrolyte and a V,O,, cathode 0.5 uA could be maintained for 120 hours. 
A rather better performance was obtained from the combination Mg anode, 
PEO,:Ca(SCN), electrolyte and TiS, cathode; 0.8 uA could be sustained 
for 1000 hours. In general, for these “mixed” systems, many of which did 
not meet our minimum performance specifications (0.5 LA, 100 hours), Mg 
anodes always outperformed Zn or Al anodes in conjunction with TiS, ca- 
thodes, whatever salt was used. This even included systems for which the 
divalent cation in the salt corresponded with the anode material. 


Table 4 shows the constant load discharge behaviour for most of the studied 
cells. Part of these studies were performed at our Laboratory and others 
were performed at the School of Chemistry of the University of Leicester 
[52]. NATO (research grant RG.86/1085) and Duracell Batteries (UK) sup- 
ported part of these studies. 


Table 4 
Cells discharge at a constant load 
Cell Max. sustained current, uÃ Discharge time, hours 
Zn | (PEO),ZnCI, | Polyaniline 1 100 
Zn | (PEO),.ZnCI, | Methylpolyaniline 1 100 
Zn | (PEO);ZnCI, |PEO, TiS. 1 30 
Zn | (PEO),NH,CI.2ZnCL, |PEO, MnO, 4 50 
Zn | (PEO),KI.2ZnCI, |MnO, 1 100 43 
Zn (PEO),ZnCL, |MnO, 100 
Zn | (PEO).Mg(SCN), | PEO, V;O,; 0.5 120 
Mg | (PEO),Ca(SCN). | PEO, TiS; 0.8 1000 
Mg | (PEO),sMg(SCN), | PEO, TIS, 0.5 400 
Mg | (PEO),«Mg(CIO,). | PEO, V;O,. 400 
Zn| (PEO),«Zn(CIO,). |PEO, TiS, 0.5 300 
Zn | (PEO),.Zn(CIO,), | PEO, V;O,, 0.5 100 
Mg | (PEO),.Ma(SCN). | PEO, TIS, 1000 
Mg | (PEO),-Mg(SCN). | PEO, TIS, 0.8 3000 
Mg | (PEO),-Mg(SCN), | PEO, TiS, >5 350(15ºC) 


7. CONCLUSION 
In recent years, investigations in the field of solid state ionics have been 


concentrated mainly on the fast ion conductors and insertion compounds and 
also on the possibility of obtaining these materials in the form of thin films. 
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We are now entering the area of solid state microionics where all devices 
will be microminiaturized and participate in the evolution of microelectronics. 
Thin film technology with a view to the fabrication of microdevices has been 
widely developed. Solid state microbatteries are produced now entirely in 
situ by molecular beam deposition, thus avoiding any eventual contamination 
of the interfaces in the microgenerator cell and producing a device built by 
integrated-circuit technology and capable of being directly integrated. R.f. 
magnetron sputtering is also being used as well as the more conventional 
chemical vapour deposition and thermal evaporation techniques. Solid state 
microbatteries are now a reality. 

A future step is the incorporation of microionic devices into integrated circu- 
its. Microbatteries are available for power sources required in integrated-cir- 
cuit architecture in order to build self-powered systems. At the minimum, in 
the present stage of development, they can be used as the holding voltage 
in many memory applications. With the integrable microbatteries available, it 
is now possible to envisage the design and construction of much memory 
and system protective architecture. 

Self-powered single devices can also be envisaged. 


In a wide range of sensor applications and in particular those in which single 
isolated events must be detected, the combination of microbatteries, with al- 
most an infinite shelf life, with a sensor constitute a promising system. 


The use of the thin film solid state battery as an energy storage device 
should be interesting, especially in association with thin film photovoltaic 
cells for solar energy conversion. 


For high energy storage the thin film battery could be envisaged if produced 
in a way analogous to superlattices with an adequate interdigitalized struc- 
ture and in an appropriate series or parallel structure. 


Microionics is now entering the forum of microdevices and very-large scale 
integration, and we should certainly witness in the coming years many inte- 
resting new applications and developments. 
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RESUMO 


Faz-se uma revisão dos funda- 
mentos da cromatografia ió- 
nica com detecção espectrofo- 
tométrica indirecta, dando 
relevo especial aos aspectos 
determinantes na sensibilidade 
da detecção. 


ABSTRACT 


The basic principles of indirect 
spectrophotometric detection in 
ion chromatography are revi- 
sed focusing the factors that 
contro! sensitivity. 


Palavras chave: lon chromatography; spectrophotometric detection; indi- 


rect spectrophotometric detection 


DETECÇÃO DE SOLUTOS EM 
CROMATOGRAFIA IÓNICA 
— DETECÇÃO ESPECTROFO- 
TOMÉTRICA INDIRECTA 


J. Costa Pessoa 
Centro de Quimica Estrutural, Instituto Superior Técnico 
Original Submetido em: 27 de Janeiro de 1993 


INTRODUÇÃO 


O processo de separação dos solutos na coluna é a questão central do mé- 
todo cromatográfico, mas a detecção das espécies separadas tem uma im- 
portância comparável. Como se referiu na parte | [1], nalguns aspectos a 
detecção em cromatografia iónica é mais problemática do que em muitos 
outros métodos cromatográficos; isto é devido a que esta é complicada pela 
quimica dos processos de separação, em maior grau do que na maior parte 
das outras técnicas cromatográficas. 


Neste texto, tal como nos anteriores [1,2], procura-se evidenciar a relação 
intima que existe, em cromatografia iónica (Cl), entre os processos separati- 
vos que decorrem na coluna cromatográfica e os sistemas de detecção. 


DETECÇÃO ESPECTROFOTOMÉTRICA INDIRECTA 


O princípio de funcionamento da detecção espectrofotométrica indirecta 
(IPD - Indirect Photometric Detection ou IPC - Ind. Phot. Chromato- 
graphy) [3,4] em Cl resulta da própria natureza particular do mecanismo 
de permuta iónica. Considerando o exemplo da cromatografia aniónica, 
em que se usa como eluente uma solução de Na*E', estando a coluna 
em equilíbrio com o eluente, todos os locais de permuta da resina estão 
ocupados com aniões E. Se se injectar uma amostra contendo um soluto 
Ct'S, os aniões S serão retidos, enquanto que os seus contra-iões Ct' 
acompanham o solvente da amostra associados a iões E', até passarem 
ao detector e esgoto. Na eluição dos iões S' há uma diminuição equiva- 
lente no número de aniões E (ver fig. 1); os princípios da electroneutrali- 
dade e da equivalência entre cargas na permuta iónica assim o exigem, 
ou seja, as concentrações totais de aniões e catiões (Na*) têm de ser 
constantes. 

isto quer dizer que a concentração de S no efluente pode ser analisada in- 
directamente pelo registo contínuo do nível do ião do eluente E. Assim, se 
os aniões que se pretendem analisar não possuem uma propriedade parti- 
cular necessária à sua detecção, por exemplo absorverem radiação de um 
comprimento de onda conveniente, pode-se explorar esta “deficiência” da 
amostra escolhendo um eluente cujo ião E' absorve radiação desse compri- 
mento de onda. Uma vez que o anião S' não absorve ou absorve em menor 
grau, será detectado pela diminuição da absorvância quando a banda con- 
tendo S' passa no detector. 
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CONCENTRAÇÃO 


Fig. 1 Eluição por permuta aniónica, evidenciando-se como a eluição de um soluto S im- 
plica uma diminuição equivalente e concertada do anião do eluente E. 


Na fig. 2 apresentam-se separações com detecção espectrofotométrica indi- 
recta. Na fig. 2(A) utiliza-se como eluente ftalato de sódio, que absorve bas- 
tante no ultravioleta e a detecção é feita a 285 nm. Pouco depois da injec- 
ção aparece uma banda positiva que é devida aos aniões ftalato deslocados 
por permuta iónica pela totalidade dos aniões da amostra; estes “não são 
retidos” e acompanham os catiões e solvente da amostra ao longo da co- 
luna. Neste caso a concentração total da amostra excede a do eluente, pelo 
que o sinal é positivo; se esta fosse menor que a concentração do eluente, 
este sinal seria negativo. Na figura 2(B) utiliza-se um eluente com quatro 
grupos carboxílicos, pelo que a sua força depende do pH e pode ser ajus- 
. tada numa gama muito apreciável. 
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Fig. 2 Exemplos de cromatogramas com detecção de aniões por IPD. 

(A) Coluna: Vydac 302 IC 4.6 a - 35 “C; eluente (2.0 ml/min): solução de ftalato ácido de 
potássio (5 mM) a pH 4.0; detecção a 285 nm. Injecção de 100 ul contendo (em ppm): 
cloreto (2.5), nitrito (2.5), brometo (2,5), nitrato (2.5), iodeto (7.5) e sulfato (7.5). Nestas 
condições experimentais o sulfato é coeluído com o pico do sistema. Os limites de de- 
tecção são da ordem de 5-30 ppb para injecções de 1 ml [5]. 

(B) Coluna: PRP X-100 (Hamilton) a 25 “C; eluente (1 ml/min): solução (0.5 mM) de piro- 
melitato a pH 4; detecção a 300 nm. Injecção de 20 ul, 0.4 mM em cada um dos iões 
(7-50 ppm). Os limites de detecção (sinal/ruído = 3) são da ordem de 25 a 50 ppb [6]. 


Algumas das questões importantes que surgem na aplicação da técnica IPD 
são: 

— como se selecciona o ião do eluente, 

— quais as propriedades espectroscópicas relevantes a avaliar, 


— qual a sensibilidade do método e de que forma está relacionado com as 
condições de eluição. 


Tal como para a detecção directa, em IPD a resposta do detector é 
medida pela diferença AA (= A,uens * Ácouo) entre os sinais devidos ao 
eluente (Asjwene = Sinal na linha de base) e o sinal quando passa a 
banda do soluto (Aqui). Se o soluto está totalmente dissociado é dada 
pelas eqs. (1) (para a detecção de aniões) e (2) (para a detecção de 
catiões) [2]: 


MA =[e, —(x/v)€,. [bes (1) 


AA 13 —(x/9) E, |DC, (10) 


em que £,. (ou £,..)e £, (ou £,.) são os coeficientes de extinção para 
o soluto S“ (ou S**) e ião do eluente E” (ou E”*), respectivamente, b é o per- 
curso óptico da célula de fluxo do detector e Cs a concentração do soluto na 
solução que passa na célula de fluxo. 


Sensibilidade na detecção 


A sensibilidade depende de vários factores e à partida seria de esperar 
que será tanto maior quanto maiores forem as diferenças entre as absor- 
tividades molares dos iões do eluente e os da amostra. Como, à priori, é 
de esperar que os iões da amostra não absorvam, ou absorvam pouco, 
concluir-se-ia que, para se obter uma sensibilidade máxima, se deveria 
escolher um comprimento de onda para o qual a absorção do eluente é 
máxima. No entanto isto pode criar problemas pois o ter-se uma absor- 
ção muito elevada devido ao eluente pode implicar não se conseguir pôr 
a zero a absorvância da linha de base e/ou ter-se um ruído muito ele- 
vado. 


sendo a concentração do eluente normalmente muito maior que a dos solu- 
tos, a passagem destes no detector é revelada por variações diminutas da 
absorvância sobre um fundo relativamente elevado. Assim, o sucesso do 
método está directamente relacionado com a possibilidade de medir 
estas pequenas variações (ou sinais) na presença do ruído associado à 
resposta do detector. 


Considerando o caso da cromatografia aniónica e que, para simplicidade de 
tratamento, a banda tem a forma de um impulso de uma onda quadrada cor- 
respondendo a uma concentração de soluto constante C., a presença deste 
soluto provoca uma diminuição na concentração do eluente do mesmo tipo 
(ver fig. 3). O sinal a medir, AA, corresponde à diferença dos sinais devidos 
ao soluto e ao eluente e é dado pela eq. (1). 

As origens do ruído neste tipo de medida não são tão facilmente identificá- 
veis como sucede na detecção condutimétrica, o qual depende bastante de 
flutuações da temperatura [3]. A absorvância não depende da temperatura 
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Fig. 3 Ilustração esquemática dos aspectos essenciais para a sensibilidade na detecção 
espectrofotométrica indirecta [4]. 

(A) Forma das bandas dos solutos e do eluente, 

(B) Resposta medida incluindo a contribuição do ruido. 


mas existirá sempre um ruído que é proporcional ao valor do sinal corres- 
pondente à linha de base [2, 3, 4]: 
Ruído=NC, €E,b (3) 


em que N é a constante da proporcionalidade. 
A relação sinal/ruído será então: 


Sinal (e, —(x/v)e E» JC; 
Ruído NE E me (4) 


Considerando o caso em que o anião do soluto não absorve para o compri- 
mento de onda seleccionado, isto é, £, =(), e não considerando os sinais, 
que não são relevantes para o que se pretende mostrar, vem: 

Sinal É: 

é E, 

Ruído NC, (5) 


sendo: 


AA =(x/y)e- Cçh (6) 


Estas expressões permitem tirar as seguintes conclusões, no que respeita à 
sensibilidade do método IPD: 


— O valor de AA aumenta com o aumento da absortividade molar do elu- 
ente (£ç), 

— O valor de AA diminui à medida que a relação x/y entre as cargas do 
ão do soluto e do eluente diminui. Isto indica que devem ser evitados 
eluentes com iões de carga superior à unidade. 


— A relação sinal'ruído e, portanto, a sensibilidade melhora à medida que 
a concentração do eluente diminui. 


Há limitações importantes no que respeita à concentração do eluente a utili- 
zar. Na verdade, tal como na detecção condutimétrica sem supressor, a 
sensibilidade aumenta com a diminuição da concentração do eluente [1]. No 
entanto a diminuição de C. (ou a diminuição da carga, y) implica maiores 
tempos de retenção e um alargamento da banda. Para evitar este efeito, 
que tem como consequência uma diminuição da sensibilidade, podem usa- 
-se colunas de menor capacidade. Corre-se então o risco de sobrecarregar 
as colunas. 


Um dos grandes méritos do método IPD é o grau de liberdade que permite 
na escolha do ião do eluente. Quase que se pode dizer que qualquer ião 
que absorva no UV ou visível se pode usar, necessitando-se apenas de to- 
mar atenção às características cromatográficas do ião, por ex. se o poder 
de eluição é excessivamente grande ou pequeno. Têm sido utilizados uma 
grande variedade de eluentes para a detecção IPD de catiões e aniões. Na 
ref. 7, págs. 122-123 e na tabela 1 estão resumidas algumas das aplicações 
e condições operacionais típicas que têm sido usadas. 


TABELA 1 - Solutos, eluentes e condições operacionais típicas para a detecção 
espectrofotométrica indirecta de alguns aniões e catiões (da ref. 2, pág. 355) 


Soluto(s) Eluente” Comprimento Limite de Ref. 
de onda detecção 
(nm) 

ANIÕES 
Aniões comuns NaNO, 247 1-10 ppm 47 
Aniões comuns Ftalato 300 0.1-0.3 ug 48 
CH,COO;, HCOO,, 
PO;”, CI Benzoato 254 0.1 ppm 49 
CIO,, Cl, SCN, SO”, 
Sd” Trimesato 317 0.5 ppm 50 
F, Cl, Br, SO,,8S,0,”, | Piromelitato 285 0.05 ppm 45 
Ct, Br, NO;, SO Citrato 290 0.2 ppm 51 
Cr, NO; NO; SOF Salicilato 254 1-2 ng 52 
Cl, Br, NO;, SO; p-Toluenosulfonato 280 50 ppb 53 
Cl, NO,, Br, NO; 2-Naftilaminosulfonato 300 2-10 ppb 53 
Cl, NO;, SO, [Cu(EDTA)P 290 1.6-11 ppb 54 
Ct, Br, NO;, SO; Tiron 290 50 ppb 55 
HPOs, Cl, NO;, SO,” Sulfobenzoato 270 300 pmol 56 
PO Po Trimesato 285 1 ppm 57 
SO", FE, CL, Br, PO; 
SO, 2,4-Dihidroxibenzoato 312 0.2 ppm 58 
SO,” 5-sulfoisoftalato 240 50 pmol 59 
CATIÕES 
lões alquilamônio 
quaternários Benziltrietilamina 270 0.6-0.8 ug 60 
Catiões monovalentes, 
aminas Benziltrietlamina 275 0.2-1 ppm 61,62 
Lit, Na*, NH; K*, Rb”, 
Cs* 4-Metilbenzilamina 262 0.2-10ppb 63 
Lit, Nat, NH.*, K*, Rb”, 
Cs 2, 6-Dimetilpiridina 269 0.4-11 ppb 63 
Lit, Na*, NH;, K*, Rb”, 
Cs* 2-Metilpiridina 262 0.4-11 ppb 63 
Lit, Na”, NHy, K, Rb”, 
Csº 2-Feniletilamina 257 0.7-14ppb 63 
Nat, K Piridina*.Cl 260 0.2 ppm 64 
Nat, NH,*, K* Cuso, 252 10 ppm 5 
Na*, K*, Rb*, Cs*, Ca”, 
Mg? CusSO, 220 30-300 ppb 65 
Na*, K*, Rb*, Cs*, Ca?, 
Mg? Ce(SO,), 254 4 ppb 66,67 
Na*, K*, Ca”, Mg* CoSO, 218 7-34 ng 68 
Ca**, Mg” [Co(en).(acac)]”* 254 10 ppb 69 
Ca”, Mg” [Cu(EDTA)|* 290 1.7-2.3 ng 54 
Mg*, Ca**, Srº*, Ba” 2-Feniletilamina 257 1-58 ppb 63 
Mg”, Ca*, Sr, Ba” 4-Metilbenzilamina 262 1-58 ppb 63 
Mg?*, Ca”, Sr”, Ba? [Cultrien)) 255 0.2-10 ng 70 


* en = etilenodiamina; acac = acetilacetonato; trien = trietilenotetraamina 
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Outra vantagem do método IPD é a versatilidade na análise de amostras 
muito concentradas em que a sensibilidade não é um factor muito relevante, 
mas interessa simplificar o método de análise. Nestes casos pode ser feita a 
análise sem diluição das amostras desde que se usem colunas de maior ca- 
pacidade e eluentes mais concentrados. 

A determinação da absorção do eluente com muito boa precisão é um 
factor determinante em IPD. A possibilidade de a absorvância devida ao 
eluente ser excessivamente elevada para que essa determinação precisa 
seja viável não deverá constituir um problema pois variando o comprimento 
de onda é normalmente possível seleccionar um valor para o qual a absor- 
vância esteja na gama 0.2-0.8, em que os erros na análise por espectrofoto- 
metria estão reduzidos ao mínimo [3,4,8,33]. 


Em termos gerais pode dizer-se que o método IPD, tendo algumas caracte- 
rísticas semelhantes ao da detecção condutimétrica sem supressor, é cerca 
de uma ordem de grandeza mais sensível [7], apresentando limites de de- 
tecção comparáveis aos da detecção condutimétrica com supressor. O facto 
de o detector também poder ser usado na detecção para cromatografia li- 
quida convencional pode também ser relevante em muitos casos. Uma das 
exigências importantes ao nivel do equipamento, é que exista a possibili- 
dade de compensar (pór a zero) os valores de absorvância elevados que 
ocorrem para o eluente. Aqui também se verifica um paralelismo com a de- 
tecção condutimétrica sem supressor, em que é frequente ter-se de com- 
pensar valores de condutividade do eluente relativamente elevados [1]. 


Aspectos particulares na detecção de aniões 


Referiu-se o paralelismo que existe em muitos aspectos entre o método IPD 
e a detecção condutimétrica sem supressor. Sucede ainda que muitos dos 
aniões que têm condutividades iónicas a diluição infinita baixas e que favo- 
recem uma boa sensibilidade na detecção condutimétrica, têm também valo- 
res elevados de coeficientes de extinção. Isto quer dizer que muitos dos elu- 
entes e condições operacionais são comuns aos dois métodos, o que não 
só facilitou o desenvolvimento do método IPD como faz com que seja possí- 
vel que os dois tipos de detecção sejam utilizados em sequência [8]. 


O pH do eluente é uma variável importante em IPD. Não só influencia a 
carga do anião deslocador como pode afectar a absortividade molar do elu- 
ente. Além disso a utilização de um pH pouco adequado pode ter como con- 
sequência uma diminuição da sensibilidade do método. Por exemplo, valo- 
res de pH excessivamente elevados, quando se usam como eluentes sais 
de ácidos aromáticos em separação de aniões, podem ter efeitos negativos 
porque os idbes OH vão competir com os do eluente e contribuir para a elui- 
ção dos aniões do soluto. Como os iões OH' não absorvem radiação, o sinal 
observado vai diminuir. É por esta razão que a sensibilidade na detecção, 
quando se usa o ftalato como eluente, começa a diminuir para pH > 7.0 
[8.34], para o qual o ftalato está completamente dissociado. 


Sucede ainda que eluentes com aniões completamente ionizados têm me- 
nos tendência a formar picos do sistema do que os que estão parcialmente 
dissociados [1,3,8,34,35]. Isto significa que a escolha do pH do eluente terá 
de obedecer a um compromisso entre o ter-se uma sensibilidade máxima e 
a eliminação dos picos do sistema. 


Aspectos particulares na detecção de catiões 
As bases aromáticas constituem uma classe importante de eluentes que têm 


aplicação no método IPD para a análise de catiões e à primeira vista pode- 
ria parecer que seriam os eluentes ideais com esta finalidade, pois consti- 


E T—O o O a me 


tuem a contrapartida em cromatografia catiónica dos ácidos aromáticos em 
cromatografia aniónica. 

Esta previsão é parcialmente verdadeira mas a aplicação destes eluentes 
deve ser feita criteriosamente. Na verdade, quando se utilizam bases aro- 
máticas a pHs relativamente baixos, por exemplo o ácido picolínico a pH 
3.5, verifica-se que a sensibilidade é muito má [36]. A origem deste facto re- 
sulta de que a eluição dos solutos catiónicos está a ser feita pela base aro- 
mática protonada e pelos iões H* presentes. Como o ião H'* não absorve ra- 
diação e a sua presença na banda do soluto implica uma diminuição da 
concentração da espécie que absorve, a sensibilidade diminui. 

Nestas condições, se se pretende uma boa sensibilidade, só as bases aro- 
máticas que estão significativamente protonadas a pH neutro têm aplicação 
no método IPD para a análise de catiões. Como exemplos de bases aromá- 
ticas adequadas podem referir-se [8,25]: 2,6-dimetilpiridina (pKa - 6.7), 2- 
-metilpiridina (pKa - 5.9), 2-feniletilamina (pKa -9.8) e a 1,4-fenilenodiamina 
(pKa - 6.2). 

Na tabela 1 e na fig. 4 ilustram-se aplicações do método IPD para a análise 
de catiões mono e divalentes utilizando como eluentes soluções de bases 
aromáticas. 


[2s=10+UA 


Mo” 


Fig. 4 Detecção espectrofotométrica indirecta de catiões utilizando bases aromáticas 

como eluentes e uma coluna Waters IC-PAK C (50 x 4.6 mm) at = 35º € [25]. 

(A) Cromatograma para catiões monovalentes; eluente: solução de 4-metilbenzilamina (0.1 
mM) a pH 6.92 e injecção de 15 ul de uma solução contendo 2x10º M de cada ião. 

(B) Cromatograma para catiões divalentes; eluente: solução de feniletilamina (10 mM) a 
pH 5.49 e injecção de 60 ul de uma amostra contendo 10º M de cada ião. 


Nos exemplos da figura 4 o detector UV está colocado em série com um de- 
tector condutimétrico sem supressor. Verificou-se [25] que de uma forma ge- 
ral na separação .de catiões monovalentes os limites de detecção são me- 
lhores com o detector condutimétrico no caso do NH,*, Rb* e Cst e 
aproximadamente os mesmos no caso do Li, Na'* e K*. No caso concreto 
das soluções de 4-metilbenzilaminas os limites de detecção são, em ppb 
[25]: Lit (0.3; 0.2), Na* (0.8; 0.8), K* (0.4; 0.7), NH, (0.6; 1.3), Rb* (1.3; 4.1) 
e Cs* (2.5; 10.2), indicando-se em primeiro lugar o valor para o detector 
condutimétrico. 

No caso dos catiões divalentes e usando a 2-feniletilamina como eluente, 
verificou-se que com o método IPD se obtiveram melhores condições de de- 
tecção do que com a detecção condutimétrica sem supressor. Mais concre- 
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tamente, os limites de detecção obtidos foram, em ppb [25]: Mg* (9; 1.3), 
Ca”? (19; 3.8), Sr” (48; 12) e Ba** (82; 58), indicando-se em primeiro lugar o 
valor para o detector condutimétrico. 

No método IPD para a análise de catiões também podem ser usados eluen- 
tes com base em soluções de sais inorgânicos (ex. Cu(ll), Co(ll), Ni(ll), 
Ce(Ill)) e complexos iónicos, ex. [Cu(en),]*. Na tabela 1 e fig. 5 apresen- 
tam-se exemplos de análises de catiões com estes eluentes. 
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Fig. 5 Exemplos de separação de catiões com detecção espectrofotométrica indirecta utili- 

zando eluentes inorgânicos. 

(A) Coluna: Interaction ION-210, eluente (1 ml/min): soluções (0.1 mM) de sulfato de 
Ce(Ill) e detecção a 245 nm. Injecção de 20 ul contendo (em ppm): Na (1.6), K (2.1), 
Rb* (7.1), Cs' (12), Mg” (2.5) e Ca” (1.6). Os limites de deteção foram (em ppb): Na* 
(4), K* (20), Rb' (100), Cs* (200), Mg” (50) e Ca” (90) [28]. 

(B) Coluna: Mitsubishi MCI CPK-08 a 30 “€, eluente (1.0 ml/min): solução (10 mM) de 
[Cu(en).]SO, e detecção a 295 nm. Injecção de 20 ul contendo (em ng): Na* (230), K* 
(390) e NH,' (180). Os limites de detecção determinados neste trabalho [16] (sinal/ru- 
ido = 2) são inferiores a 10 ng do ião na amostra injectada, ou seja, aproximadamente 
da mesma ordem de grandeza do obtido com o método correspondente a (A). 


Calibração em IPD 


É prática comum em IPD inverter a polaridade de maneira a que os croma- 
togramas apresentem o seu aspecto usual. Nas figuras apresentadas isto 
não é feito para que fique bem evidenciada a natureza do fenómeno. 

A calibração é feita da maneira usual injectando padrões de concentração 
conhecida e medindo as 'alturas' das bandas ou a sua área, A questão im- 
portante é que para muitos iões a área da banda é independente do ião in- 
jectado e depende apenas da quantidade injectada. Isto constitui uma con- 
sequência natural da natureza do método de detecção e só não é válido 
quando os solutos absorvem radiação do comprimento de onda utilizado na 
detecção. Para estes será necessária uma calibração individual [3,4]. 


Detecção espectrofotométrica indirecta e evolução da cromatografia iónica 


A detecção espectrofotométrica indirecta é uma técnica perfeitamente esta- 
belecida para a análise de aniões e em menor grau para catiões. O mesmo 
tipo de conceito tem sido aplicado à detecção condutimétrica (ex. ref. 1 e 
37), indice de refracção (ex. ref. 10) e fluorescência (ex. refs. 8, 38-42), fos- 
forescência (ex. refs 43 e 44) e absorção atômica (ex. ref. 8 e 45). 


A evolução da Cl tem sido orientada em boa parte para a determinação si- 
multânea de aniões e catiões, tendo sido abordadas principalmente as 
áreas: 


(a) 


(b) 


2 


(d 


erica 


Od 


IPD utilizando um único detector, colunas de permuta aniónica e catió- 
nica ligadas em série e procedendo a uma única injecção de amostra. 
É possível obter-se o cromatograma conjunto da resposta dos aniões e 
catiões e os cromatogramas individuais (por desconvolução do primeiro) 
quer haja coeluição destes ou não (ex. refs. 3 e 46). 

Desenvolvimento de técnicas de determinação de catiões e aniões com 
colunas em série mas em que, através de um sistema adequado de ac- 
cionamento de válvulas, o eluente só passa por uma coluna de cada 
vez (ex. ref. 47). Tal como em (a) procede-se a uma única injecção de 
amostra e apenas se usa uma bomba e um detector (condutimétrico ou 
espectrofotomêétrico). 

Desenvolvimento de colunas “mistas” para separar aniões e catiões e 
das técnicas adequadas, nomeadamente: eluentes, controlo de pH e 
força iónica para uma separação adequada (ex. refs. 48-51). 
Desenvolvimento de novas fases estacionárias quer para cromatografia 
aniónica, quer catiónica, com detecção condutimétrica ou espectrofoto- 
métrica (ex. ref 52). 


esenvolvimento de cada uma destas áreas sai no entanto fora do âmbito 


e objectivos deste texto. 
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ABSTRACT 


This text intends to be a brief 
description of the existing rela- 
ton betwen classical and 
quantum mechanies in Heisen- 
berg's early formulation. The 
historical evolution is closely 
followed although the text do- 
esn't intend to be a historical 
essay. 

The purpose of this text is also 
to present quantum mechanics 
from a different perspective 
and, perhaps, contribute for a 
better understanding of this 
fascinating subject. 


A MECÂNICA QUÂNTICA 
DE HEISENBERG 


“ a small path that led from the darkness towards the light of a new physics...” 


João Paulo Neves Monteiro dos Santos” 
Original Submetido em: 23 de Novembro de 1994 


1. INTRODUÇÃO 


O objectivo deste texto é apresentar a mecânica quântica tal como ela sur- 
giu primeiramente. Não será portanto referida a formulação posterior de 
Schrôdinger que lhe é equivalente nem o trabalho de interpretação da nova 
mecânica que se lhe seguiu. 


2. A CINEMÁTICA QUÂNTICA 


1. TEORIA DA RADIAÇÃO: MECÂNICA QUÂNTICA VS. MECÂNICA 
CLÁSSICA 


Historicamente a teoria quântica da radiação baseia-se nos três pressupos- 
tos seguintes: 


|-—- Um sistema limitado é estável apenas para um conjunto discreto de 
energias às quais correspondem estados ditos estacionários. 


| - Uma vez que cada estado estacionário possui energia constante, ape- 
nas podem ocorrer transferências de energia com o exterior como re- 
sultado da transição do sistema entre dois estados estacionários. 

IH — A radiação emitida ou absorvida em cada transição é monocromática 
sendo a sua frequência dada por £,-£, -=hw em que E, ck, são 
as energias dos estados n e n' respectivamente. 


Classicamente, dado um sistema periódico as suas variáveis podem ser es- 
critas em série de Fourier. 


+ 
a PDT 
A= > Ae 


t=—=s 


Um sistema deste género emite ou absorve radiação de frequência TO. 
Tudo se passa como se o sistema fosse substituído por um conjunto de os- 
ciladores harmónicos oscilando com frequências TX, cada um contribuindo 
com a sua frequência para o espectro. 


Podemos ver o sistema quântico como um conjunto de osciladores virtuais de 
frequências correspondentes a cada risca do espectro. Tudo se passa então 
como se uma transição entre dois níveis fizesse ressoar um oscilador virtual 
com a frequência adequada responsável pela emissão (absorção) de energia. 

Fisicamente, a origem das frequências dos osciladores virtuais é completa- 
mente diferente em mecânica clássica e em mecânica quântica. Classica- 


"* Bohr, N. On the theory of line spectra in [Van der Waerdenl 
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mente a frequência da radiação está estritamente relacionada com o período 
da trajectória. Na teoria quântica está relacionada com o espaçamento entre 


níveis de energia. Classicamente temos, usando variáveis acção ângulo 
DT = id T 
dl 


e quanticamente 


co(n + cm=[H(n+ )- H(n)| (1) 


sendo w(n+7T,n) a frequência emitida na transição do nível n+T;, para o 
n. Desenvolvendo o segundo membro desta equação em série de Taylor 
obtém-se 


| dH | dH 
o(n+ Tun)=— T+o(T)= 


h On d(nh) EA) 


A igualdade das expressões clássica e quântica no limite obriga a que 
ande ed 
on (2) 


para C=() estas são as conhecidas regras de quantização de Bohr- 
-Sommerfeld. Em rigor pode-se provar que € = *, 

É importante notar que a passagem para a mecânica clássica se faz à custa 
da transformação duma diferença numa derivada. Este processo é perfeita- 
mente geral. 

Uma diferença essencial entre sistemas clássicos e quânticos é que em me- 
cânica clássica o sistema emite continuamente radiação enquanto um sis- 
tema quântico apenas emite radiação quando se dá uma mudança do nível 
de energia. De facto, para " grande os níveis de energia estão muito próxi- 
mos e a interacção do sistema com o exterior (que existe sempre mesmo 
nos sistemas mais bem isolados), que pode ser perfeitamente desprezável 
para outros efeitos, torna-se demasiado grande em comparação com a dis- 
tância entre níveis o que origina a instabilidade prevista pela teoria clássica. 
Além da frequência há que considerar também a intensidade da radiação. 
Em mecânica quântica esta é dada por nho sendo n o n.º de fotões de 
frequência «w emitidos (absorvidos). Como em cada transição é apenas emi- 
tido um fotão, a intensidade está relacionada com o número de transições 
ou, em última análise, com a probabilidade de uma transição ocorrer. Em 
mecânica clássica a intensidade está relacionada com a amplitude de osci- 
lação do oscilador virtual associado à transição e responsável pela emissão 
(absorção) da radiação. É portanto razoável supor que esta amplitude esteja 
de alguma forma relacionada com a probabilidade de transição. 


2. AS AMPLITUDES DE TRANSIÇÃO 


Em mecânica quântica as grandezas físicas, que em mecânica clássica esta- 
vam estritamente relacionadas com as órbitas periódicas (as trajectórias) pas- 
sam a depender de dois estados estacionários. É pois razoável representar 
uma grandeza física como um conjunto de osciladores virtuais, ou seja, pelas 
amplitudes de transição Ae” dependentes de dois estados estacionários, 
em que A representa a probabilidade de transição e w a distância entre os 
níveis. Ao fazê-lo abandonamos o conceito de trajectória para nos fixarmos 
apenas no que é observado experimentalmente, as transições entre níveis. 


* Ver por exemplo [Landau] 


Representemos as amplitudes por 4, e w,= ww. em mecânica clássica e 
por A(nn-co) e w(n,n—a) em mecânica quântica. Põe-se agora o pro- 
blema de encontrar as leis a que obedecem estas quantidades de transição. 
Façamos uma analogia com o caso clássico: A soma entre duas grandezas 


escritas em série de Fourier é dada classicamente por 


à 


A+ B= SUA, + B, Jetro = SKA i. B), give 


[=—a [I==00 


pelo que quanticamente é natural escrever 
A+B-[A(n,n-0)+B(n,n- 0) =(A+B)(n,n-— oe! (3) 


Identicamente teremos, para a derivação em ordem ao tempo 


À — SA(nsn - À dad = io(n,n- o)A(n,n— ate! 


di (4) 


O problema da multiplicação é mais complicado. Consideremos duas séries 
de Fourier 


dC] E] 
A N MoTA a IEL 
A = ) Ae B= > B.e 
f=-—»e 


O seu produto dá 


AB=5 DA B ÃO = Z|Z Age as = > A, ii = 5 (AB), ee 
Tr ql clA “ 


Quanticamente deveriamos ter algo como 


AB - (AB)(n,n Ro À Fi se E fa(nn -o +B)B(n,n O ic, 
f 


No entanto esta expressão não pode estar certa pois as exponenciais são 
incompatíveis. Classicamente as exponenciais eram compatíveis porque é 
válida a relação 


+to,=o(a-B)+o.B=0, 
que permitia às exponenciais dentro do somatório não dependerem de Pp. 


A relação análoga em mecânica quântica é o princípio de combinação de 
Ritz que é uma consequência directa da estrutura de níveis: 


o(nn-B)+o(n-Bn-a)=o(nn-o) 


Usando esta relação podemos então escrever 


AB - > A(n,n Es B)B(n E B, pr ooo tet” B rem n- Ba) a (AB)(n, ni > À sda 
i (5) 


Um facto novo bastante importante surge com esta regra de multiplicação. 
A multiplicação deixa de ser comutativa, ou seja, em geral AB-BA 0. 
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3. SISTEMAS QUASIPERIÓDICOS 


Num sistema quasiperiódico com n graus de liberdade as grandezas fisicas 
podem ser escritas na forma [Goldstein] 


A relação com a mecânica quântica faz-se exactamente da mesma forma. 
Tem-se então 


LoH a qui | 
bite =) EA o(m, + tn)=[H(m,+7,) H(n,)] 


obtendo-se desta vez as regras de quantização na sua forma mais geral 

!, =hn,, Como consequência destas últimas relações passam a existir mais 

que um número quântico a caracterizar o sistema, pelo menos no limite 

clássico, o que tem uma consequência importante: se as frequências não 
k 

forem independentes, isto é, se se puder ter 3 to, = então ocorre uma 


1= 
degenerescência, isto é, o sistema quântico possui dois níveis distintos 
com a mesma energia. Tais sistemas dizem-se degenerados. 


4. CONSTANTES DO MOVIMENTO 


As constantes do movimento adquirem uma forma particularmente simples. 
Para que as amplitudes de transição A(n,n— Je“ “ sejam independen- 
tes do tempo é necessário que 


o(nn-c)=0 ou A(nn-e)=0 


Ora para sistemas não degenerados w(nn-0)=0€>0=0U pelo que ape- 
nas os termos da forma 4(1,n) são diferentes de zero. Este resultado pode 
ser estendido ao caso de sistemas degenerados-. 


3. A DINÂMICA QUÂNTICA 


Até agora temos estado preocupados com a cinemática quântica, ou seja, 
o modo como as grandezas físicas são representadas na mecânica quân- 
tica. Falta descobrir as leis físicas que governam as quantidades que estu- 
dámos. 

Consideremos agora a equação (4). Usando (1) e notando que apenas os 
H(m,n) são diferentes de O obtemos a importante relação 


Ã= (Hinn) —-H(n-a,n— e)JA(n,n— ae! = THA — AH) 
Re | 
Torna-se portanto de importância crucial a compreensão do significado de 


expressões da forma AB- BA, expressões a que chamaremos daqui em di- 
ante comutador. 


5. AS REGRAS DE QUANTIZAÇÃO DE BOHR-SOMMERFELD 
MECÂNICA QUÂNTICA 


Calculemos as variáveis acção em mecânica quântica 
1=-—6 pag =—— | púdi 
=-—— O pdq = — é 
2x pes 2x "dá 


Temos então, usando a equação (4) 


I(n+a,n)= A “SZp(n+o,n+B)lio(n+Bnag(ln+Bn)fe' ar 
27 1) ; 


Em mecânica clássica o integral é tomado sobre um periodo. Este conceito 
não pode ser facilmente transferido para mecânica quântica pelo que a ex- 
pressão acima apenas pode ser considerada no limite clássico. Aí, a relação 
2K 
o(n+B.n)= B.  , permite-nos considerar o período como sendo 1 = EE 


Integrando, todos os termos excepto os que não contêm o tempo dão zero 
devido à periodicidade da exponencial. Obtemos então 


Hn.n)= — R Sio.p. p(n.n+B)a(n+B.n)dt = DB. p(n.n+ B)g(n+ B.n) 
A » p 


Segundo a equação (2) temos então 


É (1,0) =H= De slot + BJg(n +B,n)] 


mn 


Tranformando as derivadas em diferenças obtemos finalmente 


-ih= 5 [o(n.n+B)a(n+B.n)- p(n—B.n)a(n,n - 8) =(pg-qgpXn,n) 
B 


Vejamos agora que todas as amplitudes (pg -gp)(n+0,n) se anulam: 


Para um hamiltoneano da forma H = SA V(q) é fácil de ver que 
d oV 


= (pa)= dat pã=-—>q+pÊ = GEL reapaip= Ela 
dt dg m dq m * aadda 


" 
logo a AP4 -qp)=0, pg-qp é uma constante do movimento e é válida a 
importante relação 
pg-gp=-ih (7) 


Este resultado pode ser estendido a hamiltoneanos mais gerais. 
6. O LIMITE CLÁSSICO DO COMUTADOR 


Dadas duas grandezas físicas A e B calculemos 4B-BA no caso unidi- 
mensional 


(AB-BA)nn-a-B)=A(nn-a)B(n-a,n-a-B)-B(nn-B)A(n-Bn-a-pB) 


* Trata-se do problema da diagonalização duma matriz com valores próprios de mul- 
tiplicidade superior a sem. 
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onde por simplicidade se suprimiu os somatórios. Prosseguindo obtemos 


. Ann 0)-Aln-Bn-ae-B) e (n,n—B =| B(n. n-B)-B(n—- a, n-B-c)jA(n. n—€) 


Podemos agora passar as diferenças a derivadas: 


A(n,n-0)-A(n-Bn-a-B)= BS (n-pun-a- B) 


Ficamos então com 


oB | 
GA —(n-Bn-a-B)B.B(nn-B)->—(n-on-B-ojo A(nn-—o) 
on dn 
Como as amplitudes de transição são da forma Ac“ com w=wt a variável 
angulo obtemos 


o. A(nn— 0)= nad -— 0) 
dw' 


e substituindo acima temos 


À Son — B, H — ap) SÉ (nun -B)- So = :68,h B = a) (nn es | 


H W CM CW 


Usando as regras de quantização de Bohr-Sommerteld temos 


mn Sa n-0)S(n-an-a- )-SE (nn 6) SH (n- Bina ) 


Obtemos finalmente a relação pretendida 


AB— BA = ih—— —> — —— — -=Ih,A4,B 
Po a (8) 


em que [4.B) é o parênteses de Poisson. 
Para sistemas com vários graus de liberdade 4 e B dependem de mais nú- 
meros quânticos. As contas prosseguem como anteriormente com a res- 
salva que agora se tem 
k 
Aln,n,— 6, Mi Ph , hai 
(n0m-0,)-A(m Bom 0-8) DB So E 


de tal forma que no final se obtem 


E E A,B| 
= ow, ol, dw, dl, 


AB — BA - ns 


O limite clássico do comutador é portanto o parênteses de Poisson. 
7. O SIGNIFICADO DO COMUTADOR 
A relação entre o comutador e o parênteses de Poisson vai muito além do 


visto na secção anterior. De facto, classicamente temos A=[A,.H), expres- 
são em tudo análoga à equação (6). Por outro lado a expressão (7) tem o 


seu equivalente clássico na igualdade (qo D;|= 0, O parênteses de 
Poisson não é portanto apenas o limite clássico do comutador. O comutador 
é o equivalente quântico do parênteses de Poisson. Assim, a dinâmica 
quântica é em tudo semelhante à clássica. No entanto, a representação das 
grandezas, a cinemática, é algo completamente novo: descreve apenas 
aquilo que é observado experimentalmente — transições entre níveis — 
tendo desaparecido o conceito base da mecânica clássica, a trajectória. 


4. CONCLUSÃO 


Foi assim que surgiu primeiramente a mecânica quântica. Seguiu-se um 
longo trabalho de compreensão da teoria e das diferenças profundas que a 
separam da mecânica classica do qual podemos destacar a interpretação 
probabiliística da função de onda, de Born, as relações de incerteza de Hei- 
senberg e, mais recentemente, a desigualdade de Bell. 
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RESUMO 


Apresenta-se a equação de 
Duffing e o método de Liapu- 
nov-Schmidt. Aplicando este 
último mostra-se a existência 
de pontos homoclínicos trans- 
versos na equação de Duffing 
perturbada e consequente- 
mente a existência de compor- 
tamento caótico. 


EQUAÇÃO DE DUFFING 
E METODO DE 
LIAPUNOV-SCHMIDT 


PONTOS HOMOCLÍNICOS TRANSVERSOS E CAOS 


Diogo Luís Aguiar Gomes” 
Original Submetido em 7 de Dezembro de 1994 


1. A EQUAÇÃO 


A equação de Duffing surgiu originalmente no estudo de fenómenos eléctri- 
cos em 1918 (artigo original Duffing, G “Erzwungene Schwingungen bei Ve- 
rândelicher Eigenfrequenz und ihre Technische Bedeutung”, Braunschweig. 
Samm Vieweg, 1918). Ela aparece ligada a problemas de sub-harmónicas 
em linhas de alta tensão ou circuitos oscilantes com saturação magnética 
[2]. É ainda utilizada para descrever diversos fenómenos físicos como por 
exemplo a trajectória de uma partícula de massa unitária num potencial da 
forma 


(1) V(x)=-> + 


cujo hamiltoneano é H(x,X)=5-5 ++ [9]. 
Quando é introduzido atrito e um forçamento periódico, a equação de Duf- 
fing toma a forma 


(2) ps 


|5 =x-x'+A,y+A, cos(t). 


Diversas propriedades interessantes têm sido observada em sistemas que 
obedecem a esta equação. Por exemplo, variando os valores dos parâme- 
tros aparecem órbitas de períodos cada vez mais longos, sempre sub-múlti- 
plos do forçamento, até que a uma dada altura se atinge uma região de 
aperiodicidade [2]. Neste artigo procuramos apresentar uma nova prova da 
existência de comportamento caótico [4] por parte das soluções do pro- 
blema utilizando o método de Liapunov-Schmidt. Este corresponde de certo 
modo a um teorema da função implícita generalizado, que inclui como caso 
particular o método de Melnikov vulgarmente utilizado para mostrar a exis- 
tência de caos (por exemplo [3] e [4] utilizam o método de Melnikov para 
demonstrar a existência de comportamentos caóticos na equação de Duf- 
fing). 

O sistema não perturbado (4, = A, =0) possui duas órbitas homoclínicas 
da forma [4] 


(3) ole tv2sech(r) 
Polt)| | Fv2sech(t) tanh(t) |. 


* Monitor no Dep. de Matemática/IST 
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Na última secção demonstraremos que é através da perturbação destas ór- 
bitas que aparece o comportamento caótico. Assim quando este potencial é 
substituído por um da forma V(x)=--+>, com um único mínimo não apa- 
rece a bifurcação homoclínica (uma vez não existe uma órbita homoclínica 
no sistema original) não havendo portanto caos. 


2. O MÉTODO 
O método de Liapunov-Schmidt é uma das ferramentas mais poderosas 
para encontrar e descrever as soluções de problemas da forma 


(4) Ax= N(x,A) 


onde A é um operador linear contínuo de X>Z (Xe Z espaços de Ba- 
nach), N(x.A) é a parte não linear do problema (N(0,0)=0 e N,(0,0)=0) 
e À é um parâmetro (A € A, A é um espaço arbitrário de parâmetros). 


O primeiro caso que vamos analisar é quando A é não singular. Nesta situa- 
ção o método de Liapunov-Schmidt reduz-se ao teorema da função implícita 
e garante que para À numa vizinhança suficientemente pequena da origem 
existe uma única função x (4) que satisfaz (4). 

O caso seguinte a estudar é quando A é singular e existem projecções con- 
tinuas U e E tais que 


UX=X,=N(A) EZ=Z,=R(A), 


em que N(A) e R(A) são respectivamente o núcleo e contradomínio de A. 
Se dim(N(A)) - dim( Rº (4)), finitas, podemos transformar a equação (4) no 
seguinte sistema 

EAx = EN(x,À) 


(5) | ? 
(U-EJAx=(1-E)N(x,Ã) 


Se definirmos x=y+z onde yceX, e z€EX,, e atendermos à definição 
de E (notando que (/- E)JA=0) temos 


z=EN(y+z,À 
(6) | pi 


0=(I-E)N(y+2,1) 


Da primeira equação do sistema é possível encontrar uma função 


: (»,A)e X,,, única, uma vez que as condições do teorema da função im- 
plícita são verificadas. De facto existe um operador K:Z, 9X, (K=A!, 


no domínio conveniente) tal que a equação z=KEN(y+<2,A) é equivalente 
à primeira equação do sistema. Substituindo z na segunda equação do sis- 


tema, obtem-se a função de bifurcação 


(7) G(y.A)=(1-E)N(y+2"(y.A),4) 


cujos zeros (se existirem) determinam as soluções de (4). 


Temos portanto um método que permite dividir um problema de dimensão 
elevada ou até infinita num problema de dimensão menor (eventualmente de 
dimensão 1) e noutro cuja solução é garantida pelo teorema da função im- 
plicita. 


3. PONTOS HOMOCLÍNICOS TRANSVERSOS 


Definimos a aplicação de Poincaré (P) para a equação (2) como sendo a 
aplicação que ao ponto p associa o ponto P(p)=P(p.27) onde D(a.t) é 
a aplicação fluxo induzida pela equação diferencial'. A aplicação de Poin- 
caré dá-nos assim a solução da equação diferencial com condições iniciais 
p após um período de forçamento. 

Chamamos a um ponto q ponto homoclínico a p (p ponto fixo hiperbólico e 
diferente de q) se a sua órbita através da aplicação de Poincaré se aproxi- 
mar de p qualquer que seja o sentido do tempo (i.e..lim,,..P'(g)=p) é 
então imediato que q pertence à intersecção das variedades estável e instá- 
vel do ponto p. 


Se a intersecção em q das variedades estável e instável for transversal en- 
tão o ponto q será chamado de homoclínico transverso. A existência de um 
destes pontos implica a existência de uma infinidade destes (basta pensar 
nas iteradas de q). Implica também que as variedades estáveis e instáveis 
do ponto p se acumulem sobre si próprias de um modo complicado como 
mostra o seguinte lema [5] 


Lema 1 (lema da inclinação ou À-lema). Seja f um difeomortismo C' em 
R" com um ponto fixo hiperbólico p. Sejam s e u a dimensão das varieda- 
des estável (W'(p)) e instável (W'“(p)) (s+u=n), Seja então D uma bola 
de dimensão u em W'(p) e K uma bola de dimensão u que intersecta 
W'(p) transversalmente num ponto q. Então U,o f"(K) contém bolas de 


dimensão u arbitrariamente C' próximas de D. 


No caso de estarmos em presença de uma intersecção homoclínica trans- 
versa podemos escolher K de modo a que esteja em Wº(p), tal implica 
imediatamente que as duas variedades se dobram e acumulam de uma 
forma complicada como se pode ver na figura seguinte. 


De facto é possível mostrar que a aplicação de Poincaré tem um conjunto 
numerável de órbitas periódicas com períodos arbitrariamente grandes, bem 


1 


considerando a equação i=| com resS!, tal como é feito em [4] 
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como um conjunto não numerável de órbitas não periódicas e uma órbita 
densa. Por outro lado, é também possível demonstrar que, devido à trans- 
versalidade, qualquer difeomorfismo c'! próximo da aplicação de Poincaré 
possui uma estrutura de órbitas topológicamente semelhante a esta [1]. 


4. CAOS 


Da primeira parte sabemos que a equação (2) sem perturbações tem duas 
órbitas homoclínicas. Consideremos, por exemplo, a seguinte órbita 


á Pot) | v2sech(rt) 
(8) Palt) -"2sech(t) tanh(t) | 


Nesta secção vamos provar que para determinados valores do parâmetro 
esta órbita dá origem ao aparecimento de pontos homoclínicos transversos. 


Consideremos então o seguinte sistema de coordenadas, que acompanha a 
órbita homoclínica, 


ã 10) [o polt+ O) + ut) 
já vo Lu(t+ o) +) |, 
Então podemos reescrever a equação (2) na forma (4) 
(10) Ah=N(h,Ã,Ã,). 


onde h=[u v] com 


e 
| Ú 
WA ço daj= Pe -w +A, polt)+A,v+ As cos(t— R 
ou seja 
a 0 ER 
Pista PRA À, cos(t — o|* iu do 


Note-se que a fase à foi introduzida no coseno por mudança da origem dos 
tempos, simplificando ligeiramente a expressão. 

Estamos agora em condições de aplicar o método de Liapunov-Schmidt. De 
facto, se escolhermos para dominio do operador A o conjunto X das funções 
Rs Rº diferenciáveis e limitadas (com a norma |lhl| = suplh| + suplh|) e con- 
junto de chegada Z das funções R-—»Rº contínuas e limitadas (com a 
norma [q] =suplo|), A é linear e contínuo, Considerando a equação homo- 


génea 4h, =) é fácil mostrar que o núcleo do operador A é 


(11) SA) = span Ieçod | 


De facto temos apenas de verificar não existem mais soluções linearmente 
independentes da equação homogénea que se encontrem em X; como se 
trata de uma equação do segundo grau existem duas soluções linearmente 
independentes p(t) e g(t). Consideremos então a matriz fundamental 


| 2 f , t 
vo 
Poli) ql). 
Sabemos pela fórmula de Jacobi que 


det W(t) = det W(0) Di 


em que B(t) verifica ly(t)= B(t)hy(t) ou seja 


To | 
B(1) |, —-3po(t) |] 


Como o traço de B(t) é nulo, podemos escolher um factor multiplicativo 
conveniente de modo a que detW(0)=detW(r)=1, tal implica 
Poloalo)- a(l)p(t)=1. Como lim,.. Pol)=lim,.. Po(t)=0 é imediato 
que q(t) e q(t) não são limitadas não pertencendo portanto a X. 
Analisando o operador adjunto (A) e a equação homogénea 


A'v=0-Bwm=o, 


é imediato que a sua solução fundamental Q(t) é (não esquecendo que 
detW(0)=1) 


o(n)= ro prod 69 


at) Po(t) |. 


Assim, utilizando um argumento semelhante ao anterior, concluímos que 


(A) span) Pri | R(A) | 


Polt) 


uma vez que o núcleo do adjunto é o sub-espaço ortogonal à imagem do 
operador. 

Neste momento a aplicação do teorema de Liapunov-Schmidt é imediata. . 
Defina-se então a seguinte projecção sobre R(A) 


E toe 


: m te 
(12) Ez=z- AM m.z dt 


E Pol t) +oa , 
onde m(t)= a O dlP dt. 
(1) fra b | Info 


Podemos então demonstrar o seguinte teorema 


Teorema 1. 4 equação (10) tem solução única na vizinhança de = 
A=(A,.4,)=(0,0) se e só se se verificar 
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(13) Po pol Adolt)+ À cost 6))dr=0 


Demonstração: O método de Liapunov-Schmidt permite-nos transformar a 
equação inicial no seguinte sistema 


EAh = EN(h,À) 


(14) (I-E)JAh=(I-E)N(hA). 


Escrevendo h como h=y+z onde ve N(A) e ze N“(A) a função de bi- 
furcação é 


(15) G(v.0.2)= [| poltÃA, Pol!) + Às cos(t— 6) dr + O(|y|) 


Definindo 
(o. À) = l. PolO(A, do (1)+4, cos(1t IN 6))dr É 


G(6.1) =0 implica a existência de uma solução da equação (10). Desen- 
volvendo a equação anterior com auxílio de (3) obtemos a seguinte condi- 
ção 


24, |. sech(t)tanh(t) cos(t — óJdr 24, | sech?(r)tanh?(r)dt =0 


Integrando e simplificando tem-se — 42, +27 A,sech(£)sin(9)=0. Assim, 
só existe solução para 
+ 4cosh(5 
à, Acosh(£) 


À, 327 


De facto como se verifica Ma *() o teorema da função implícita garante a 
existência de um valor de é que é solução da equação G(v.9,A)=0 


qualquer que seja 4 e y numa vizinhança da origem E 


Finalmente chegamos ao teorema que prova a existência de caos para a 
equação de Duffing 


Teorema 2. A condição 5: *0 implica a existência de um ponto homocli- 
nico transverso. 


Demonstração: Segue do facto de se poder aproximar as variedades estável 
e instável da origem por 


x(t) al Polt+ 6) apo(t) 2 
(16) po j ira rof | 


com validade uniforme para a variedade instável para | € |= oo,7] e para a 
estável com te[T.+<| com T arbitrário. Assim é fácil de ver que para À, 
e 4 suficientemente pequenos a intersecção é transversal e portanto não 
é quebrada (ou seja numa vizinhança suficientemente pequena existe uma e 
uma só intersecção) por termos de 2º ordem. Tal implica a existência de um 


O TTÕo a TT a O e 


(e consequentemente de uma infinidade contaável) ponto homoclínico trans- 
verso. De facto se a intersecção não fosse transversal seria possivel pertur- 
bar o campo vectorial de modo que a função de bifurcação correspondente 
não se anulasse não existindo ponto homoclínico ou de modo a que existis- 
sem vários zeros da função de bifurcação arbitrariamente próximos; por ou- 
tro lado a condição % » () implica que qualquer que seja a perturbação su- 
ficientemente pequena continue a existir uma e uma só solução da equação 
(10) o que contradiz a afirmação anterior 


Com este último teorema demonstramos, por via analítica a existência de 
comportamento caótico para a equação de Duffing. 


AGRADECIMENTOS 


O autor agradece ao prof. Carlos Rocha as sugestões e o cuidado com que 
analisou o manuscrito. 


BIBLIOGRAFIA 


1. Chow, S. e Hale, K. “Methods of Bifurcation Theory”, Springer-Verlag 1982. 

2. Faro, A. “ Ordem e caos associados à equação de Duffing uma aproximação de 
natureza fisica”, Memórias da Academia das Ciências de Lisboa, Tomo XXVIII 
1987 

3. Greenspan, B. e Holmes, P. “Non Linear Dynamics and Turbulence”, Pitman Ad- 
vanced Publishing Program 1983 (G.l. Barenblatt and G. loss eds.) 

4. Guckenheimer, J. e Holmes, P. “Nonlinear Oscilations, Dynamical Systems, and 
Bifurcation of Vector Fields”, Springer Verlag 1983. 

5. Melo e Palis “Introdução aos sistemas dinâmicos”, Projecto Euclides IMPA 1978. 


71 


N > j > E | Ê E 
| | 
BA URSOS INFORMÁTICA Desa 
INTRODUÇÃO AO CAD INÍCIO: Abr/Maio/Set DURAÇÃO: 20H LOCAL: Guimarães, Porto, Leiria, Lisboa 
AUTOCAD Versão 13. a EE ea Ee Ê 
CURSO EM TEC. DE ENGº. ASSISTIDA POR COMP. INÍCIO: 8 Maio DURAÇÃO: 96H LOCAL: Porto 
ÁREA  [ANSYSPC Du o ss 
DE CURSO MULTIMEDIA INÍCIO: 22 Maio | DURAÇÃO:90H — LOCAL: Porto e Lisboa 
=“ 4n4 | AUTHORWARI 
CAD/CAM DESENHO MODULAÇÃO E ANIMAÇÃO 3D nícIO:22 Maio | DURAÇÃO: [20H LOCAL: Porto 
| 3DS TUDIO Versão | | 
Horário | CAD/CAM AVANÇADO INÍCIO: | Setembro DURAÇÃO: 160H | LOCAL: Porto 
Pós-Laboral AUTOCAD Versão 13 
| CURSO MODULAR APLICATIVOS AUTODESK iníCIO: 11 Setembro DURAÇÃO: 20H LOCAL: Porto 
ADE. AUTOVISION, AUTOSURF, DESIGNER 


Caso pretenda obter informações mais detalhadas acerca de alguns destes cursos, preencha este destacável e envic-o paro o nosso apartado no AT 
Porto: APARTADO 1145 - 4102 - PORTO CODEX - Receberá em sua casa ou na sua empresa os respectivos folhetos. - 


NOME: 
EMPRESA: === 1 = = —— 
MORADA: =" io a 
Erosr Lit). cs “TELEF: FAX: 


Coloque a sua movcada particular ou da empresa conforme pretenda receber a resposta em casa ou na empresa. Ministramos formação personalizada à empresas. 


CENTRO DE FORMAÇÃO PROFISSIONAL DE INFORMÁTICA - CESAI 
PORTO: Rua Ciriaco Cardoso, 186 - 4150 PORTO Telet.: 610 16 76 Fax: 610 205] DELEGAÇÃOES: Aveiro, Coimbra, Guimarães, Leiria, Lisboa, ViseL 


Somos Centro de Formação Autorizado das Seguintes Marcas: 


; ay 
ns : ai : | 
q IN Autodesk Microsoft e Sia 
Es | E E “SOLUTION PI E V] To! R - Ri Ta DT Bia 
drpatio Ehadot [ralninog Center PRE mi isisadshedis i 


FEIRA DE ESTÁGIOS E EMPREGO 


2 a 5 de Maio 1995 


JUVENTUDE 


LABELEC 


Grupo EDP 


ad 


E 


“ QUALIDADE PORTUGUESA 


PRESTÍGIO INTERNACIONAL 


a a 


e, 
SE bado et 


| 


AL. 


, 
| 


Diagnostico do estado 
de transformadores, 
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cromatografica dos 
gases dissolvidos 

no óleo. 


saio de tipo de uma 
une de isoladores 
400 RV 


Desenvolvendo actividades técnico-cientificas para o Grupo EDP e outras empresas, 
os Laboratórios LABELEC estão equipados para oferecer uma ampla variedade de 
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A principal vocação da ABB 
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Para a ABB, isto 
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Esto não é ficção cientifica. E parte do dia-a-dia 
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